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SECTION  1 


■  :  INTRODUCTION 

j 

The  objective  of  the  present  work  is  to  assess  the  feasibility  of 
studying  the  blast,  fire,  and  wind  effects  due  to  a  nuclear  explosion  in  an 
urban  setting  by  employing  scale-modeling  techniques.  Our  concentration  will 
be  on  fire-scaling.  The  blast  effects  will  be  considered  only  to  the  extent 
of  concern  to  describe  the  initial  combustibles  distribution  r  ground-zero. 
In  a  similar  way,  the  winds  produced  directly  by  the  explosion  will  be 
considered  only  as  to  their  impact  on  the  initial  fires.  As  the  blast  wave 
passes  and  the  associated  winds  subside,  the  fire  would  develop  in  a 
self-regulated  manner  involving  complicated  heat  transfer,  turbulent  buoyant 
fluid  mechanics  and  combustion. ^It  is  this  development  upon  which  we  focus 
attention. 

Existing  literature  indicates  that  knowledge  about  the  behavior  of  fire 
following  a  nuclear  explosion  is  less  complete  than  about  the  blast  and 
radioactive  fall  out.  A  nuclear  explosion  triggers  fires  in  an  urban  area  by 
two  different  mechanisms:  First,  primary  fires  are  initiated  due  to  thermal 
radiation  emitted  by  the  fire-ball  over  an  area  in  close  enough  proximity  to 
ground-zero  where  the  irradiance  is  sufficiently  high  to  heat  the  combustibles 
to  the  state  of  ignition.  These  nascent  fires  are  then  subjected,  within 
seconds,  to  the  blast  effects  of  wind,  pressure-rise  and  mechanical 
disruption.  The  consequence  may  be  either  extinguishment  or  augmentation  of 
the  fires.  Second,  secondary  fires  are  initiated  as  the  blast  damages  the 
mechanical  systems  in  and  outside  the  structures  physically.  Broken  gas 
pipes,  fallen  electrical  wires,  ruptured  chemical  storage/transportation 
equipment  and  others  constitute  the  sources  of  secondary  fires. 

Fires  associated  with  nuclear  blasts  are  characteristical ly  expected  to 
continue  burning  long  after  the  blast  effects  subside.  Once  the  primary 
ignitions  are  modified  by  the  blast  wave  and  the  secondary  ignitions  are 
provoked,  the  individual  fires  are  expected  to  grow  to  coalesce  together.  In 
this  coalescing  process,  such  phenomena  as  flashover  in  enclosures  (i.e., 
rooms),  penetration  through  walls,  fire-jumping  across  breaks  in  the  fuel  bed 
are  generally  involved  in  addition  to  fire  propagation  over  contiguous  fuel 
beds.  Once  coalesced,  the  resultant  area-fire  will  conflagrate  to  result  in  a 
mass  fire. 


SECTION  2 


NATURE  OF  THE  EXPECTED  FIRE 

Thermally  thin  combustibles  are  generally  expected  to  be  the  first 
Ignited  Items  no  matter  whether  they  are  situated  indoors  or  outdoors.  These 
Include  window  coverings,  upholstery  linings,  carpets,  papers,  fabrics, 
grasses,  leaves,  dry  shrubs,  etc.  The  criterion  by  which  thermal  thickness  Is 
judged  is  based  on  whether  or  not  internal  temperature  gradients  predominate 
within  the  solid.  The  nondimensional  parameters  addressing  this  issue  is 
known  as  Biot  number  which  Is  discussed  in  detail  in  a  later  chapter  of  this 
report.  Suffice  it  to  say,  thermally  thin  bodies  are  characterized  by  a  Biot 
number  much  smaller  than  unity. 

As  energy  Impinges  on  the  solid  surface,  part  of  it  is  reflected,  part  is 
lost  by  convection  and  the  rest  Is  employed  to  heat  the  body  by  conduction  to 
ignition.  The  heated  body  may  also  suffer  reradiative  loss.  If  the  intensity 
of  irradiance  is  high,  heating  to  ignition  occurs  so  fast  as  not  to  provide 
time  to  develop  the  natural  convective  heat  loss  process.  The  ignition 
process  has  been  studied  extensively  by  experiment  (Martin  [11)  and  by 
analysis  (Kanury  [2],  fiandhi  [3]). 

Once  the  fire  is  thus  initiated,  it  will  grow.  If  the  initial 
ignition  is  indoors,  the  fire  growth  process  is  influenced  by  the  room 
configuration.  This  influence  is  weak  initially  and  strong  in  later  stages 
leading  towards  what  is  known  as  flashover.  Room  fire  flashover  is 
extensively  studied  experimentally,  analytically  and  by  scaling  (Kanury  [41). 
an  interesting  additional  feature  of  flashover  in  rooms  exposed  to  the 
fire-ball  is  the  role  played  by  the  radiant  deposition  of  energy  on  the 
unignited  combustibles  of  the  room.  This  role  is  surmised  by  Martin  to 
expedite  the  arrival  of  flashover.  The  scaling  considerations  of  flashover 
will  be  discussed  later  in  detail. 

If  the  fire  were  in  the  open,  such  as  in  a  grassy  field  or  in  a  forest, 
the  initial  fire  growth  follows  the  pattern  of  a  conflagration  at  a  small 
scale.  The  scaling  rules  governing  this  growth  are  also  discussed  later  from 
the  basis  of  the  existing  knowledge  on  fire  spread  and  forest  fires. 

As  the  fire  from  the  initial  ignitions  thus  attempts  to  evolve,  the  blast 
would  soon  arrive  to  disrupt  it.  The  effect  of  blast  impingement  is  known  to 
be  towards  extinguishing  the  gas-phase  flames  and  augmenting  glowing-type 
solid-phase  (or  near  surface)  phenomena  (Kanury  [5]).  Even  more  importantly, 
the  blast  wave  would  disrupt  the  burning  fuel  elements  perhaps  to  a 
geometrical  configuration  more  conducive  to  intensified  burning  and  would 


produce  secondary  ignitions  to  augment  the  fire. 

Finally,  as  the  fire  flashes  over  in  rooms  and  penetrates  through  walls, 
fire  spread  from  structure  to  structure  becomes  relevant  in  the  development  of 
the  conflagration  and  mass  fire.  Thermal  radiative  heat  transfer  becomes 
crucial  In  this  spread  process.  Scaling  considerations  of  both  the  radiative 
(conflagratory)  fire  spread  and  the  mass  fire  will  be  subject  to  discussion  in 
later  sections. 
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DIMENSIONAL  REASONING 

As  one  embarks  upon  the  task  of  developing  a  quantitative  understandi ng 
of  the  process  of  inception  and  growth  of  an  urban  fire,  one  encounters  a 
number  of  choices  as  to  the  level  and  detail  of  study.  Since  the  fire 
phenomenon  Is  a  complicated  mosaic  of  several  component  processes,  one  can 
choose  an  isolated  process  for  a  detailed  theoretical  and  experimental  study. 
The  more  the  numbers  of  interacting  component  processes  considered,  the  less 
aminable  the  study  becomes  to  a  detailed  analysis.  But  because  one  is  often 
concerned  with  the  behavior  of  a  real  fire,  albeit  complex,  one  is  faced  with 
the  need  to  synthesize  a  number  of  nonlinear  processes  into  a  comprehensive 
whole.  Modern  high  speed  computers  have  become  extremely  useful  in  this 
synthesis.  Experimental  methods,  however,  leave  much  to  be  desired  in  their 
sophistication,  for  real  fires  are  'dirty'  and  perhaps  will  always  remain 
benefited  from  the  modern  developments  in  instrumentation  technology.  What  is 
possible,  however,  is  powerful  data  acquisition,  quick  analysis  and 
instructive  presentation. 

While  endeavors  in  theoretical  analysis  of  fires  will  become  gradually 
more  sophisticated,  experimentation  will  continue  to  provide  not  only  the 
input  data  required  in  theories  but  also  the  output  to  validate  theories. 

The  terms  'modeling'  (broadly)  and  'scaling*  (specifically)  are 
traditionally  used  to  describe  the  endeavors  associated  with  attempts  to 
design  experiments  on  a  systematic  foundation  and  to  present  the  data  from  a 
limited  number  of  experiments  in  a  general i zable  form  suggestive  of  the 
involved  mechanisms.  Modeling  (scaling)  is  a  valuable  technique  to  understand 
the  outcome  of  the  prototype  experiment  which  is  generally  prohibitive  in 
cost,  material,  and  time,  or  when  it  involves  scales  of  time,  mass  and  size 
that  are  difficult  to  be  contrived  in  an  experimental  set-up.  Simulation  and 
analog  modeling  techniques  are  capable  of  yielding  a  convenient,  controlled, 
inexpensive,  and  meaningful  model  experiment  in  which  the  phenomena  of 
interest  can  be  highlighted. 

Note  that  the  term  scale  model  implies  more  than  the  general  connotation 
of  miniature  hardware.  The  des.gn  of  scale  model  experiments  requires  an 
Identification  of  the  constraint  parameters  which  have  to  be  preserved 
invariant  betweent  the  model  and  prototype  experiments.  These  parameters  are 
nondimensional  ratios  of  fluxes,  forces,  lengths  or  of  time  constants. 
Identification  of  these  parameters  generally  requires  an  enumeration  of  the 
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involved  variables  but  not  a  quantitative  relation  among  them. 

Scale  model  studies  of  any  physico-chemical  system  offer  one  or  more  of 
the  following  four  advantages:  (a)  attack  is  possible  on  such  new,  unexplored 
and  comple:'  problems  which  are  not  amenable  to  an  analytic  formulation  and 
solution;  (b)  transformation  to  physically  manageable  geometric  proportions  is 
possible  of  such  systems  which  are  too  extensive  to  permit  direct 
experimentation;  (c)  time  scales  can  be  expanded  or  contracted  to  speed-up  (or 
slow-down)  the  processes  which  are  too  slow  (or  too  fast)  for  direct 
experimentation;  and  (d)  since  model  experiments  can  be  conducted  redundantly 
under  well -controlled  conditions,  the  phenomenon  under  scrutiny  can  be 
understood  much  more  deeply. 

Design  of  scale  model  experiments  requires  a  preliminary, 
order-of-magni tude  analysis  of  the  underlying  physico-chemical  mechanism(s) 
of  a  phenomenon.  As  such,  the  endeavor  of  designing  a  scale  model  experiment 
squarely  falls  on  thp  boundary  between  theoretical  analysis  and  experimental 
study. 

The  idea  of  order-of-magni tude  analyses  of  a  process  automatically 
invokes  that  the  various  terms  In  an  equation  governing  the  process  behavior 
must  all  be  of  the  same  magnitude  under  tne  extreme  conditions  of  the  process. 
This  would  imply  that  prior  to  being  able  to  conduct  an  order-of-magni tude 
analysis,  one  has  to  be  able  to  postulate  a  hypothesis  upon  the  basis  of  which 
the  concerned  process  occurs  and  to  put  this  hypothesis  in  the  form  of  a 
simple  (or  complex)  governing  equation.  A  correct  governing  equation  has  to 
satisfy  two  basic  requirements:  first,  all  of  its  terms  must  be  of  equal 
magnitudes.  This  requirement  results  from  retaining  only  those  terms  which 
represent  importance  in  the  governance  of  the  process  while  ignoring  terms 
standing  for  unimportant  or  secondary  aspects.  Recognition  of  which  terms  to 
drop  out  and  which  to  keep  is  an  important  skill  of  a  knowledgeable  physicist. 
Second,  the  governing  equation  has  to  be  dimensionally  homogeneous;  that  his, 
all  of  its  terms  have  the  same  units. 

These  two  requirements  lead  to  the  following  two  axioms  which  form  the 
foundation  piers  of  the  art  of  ‘dimensional  analysis'  (Murphy  [6]). 

1.  Absolute  numerical  equality  of  quantities  may  exist  only  when  the  two 
quantities  are  qualitatively  similar.  (That  his,  an  apple  may  be  equal  to 
another  apple  but  never  to  an  orange!) 

2.  The  ratio  of  magnitudes  of  two  like  quantities  is  independent  of  the 
units  used  in  their  measurement,  provided  that  the  same  units  are  used  for 
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evaluating  each.  (For  example,  the  length  and  width  of  a  table  are  like 
quantities;  both  have  units  of  L;  the  ratio  of  length  to  width  is  same  no 
matter  whether  one  measured  them  in  Inches,  feet,  yards,  meters  or  miles.) 

Based  on  these  two  axioms,  the  powerful  tool  of  dimensional  analysis 
evolved  to  enable  deduction  of  qualitative  relationships  among  the  process 
variables.  These  relationships  guide  the  investigator  to  design  experiments 
which  yield  maximum  generalizable  data  with  an  economy  of  time  and  resources. 
The  experimental  data  in  the  framework  of  the  qualitative  relationships  lead 
to  quantitative  results  and  even  predictive  rules. 

Thus,  to  be  able  to  make  dimensional  and/or  order-of-magni tude  analyses 
of  a  given  process,  one  first  needs  enough  of  an  understanding  of  the  process 
to  identify  the  component  forces  and  fluxes.  Any  attempt  to  make  an 
order-of-magni tude  estimate  is  invariably  a  prelude  to  deducing  the 
nondimensi onal  scale  parameters  of  the  problem  and  their  expected  quantitative 
range.  As  Segel  [7]  enumerates,  there  exist  at  least  six  different  ways  to 
obtain  the  preliminary  understanding  of  any  given  process:  (a)  utilize 
experimental  and  observational  evidence  related  to  the  process;  (b)  use 
experience  from  related  or  analogous  problems;  (c)  solve  simplified  versions 
of  the  problem;  (d)  ignore  terms  which  are  small  in  order-of-magni tude  and 
solve  the  resulting  simpler  problem;  compare  the  solution  with  experiments  and 
justify  that  the  ignored  terms  are  indeed  small;  (e)  assume  certain  scaling; 
solve  the  problem;  check  if  the  dimensionless  terms  are  of  magnitude  unity; 
and  if  any  of  them  are  not,  refine  the  assumed  scaling;  and  (f)  examine  the 
complete  (or  comprehensive)  numerical  solutions  to  identify  the  scaling 
parameters . 


SECTION  4 


CERTAIN  LIMITATIONS  IN  DIMENSIONAL  REASONING 

APPLIED  TO  THE  PHENOMENON  OF  FIRE  : 

! 

( 

The  concepts  of  dimensional  reasoning  discussed  in  the  preceding  chapter  1 

are  Indeed  powerful  in  solving  numerous  problems  in  science  and  engineering.  ! 

Much  progress  has  been  made  in  fire  studies  also.  There  exist,  however,  : 

several  limitations  of  the  modeling  technique  when  dealing  with  fire  ; 

phenomena.  It  is  worthwhile  to  formally  identify  and  constantly  remember  | 

these  limitations. 

Combustion  phenomena  involve  a  host  of  physico-chemical  processes 
including  heat  and  species  diffusion,  viscosity  effects,  buoyancy,  surface 
tension,  compressibility,  kinematics  and  dynamics,  heterogeneous  and 
homogeneous  reactions,  equilibria,  heat  and  mass  transfer  to  (or  from)  the 
boundaries,  thermal  radiation,  condensation  and  evaporation,  particle 
agglomoration,  and  others.  If  the  characteristic  rates  (or  times)  for  these 
various  processes  can  be  identified,  then  one  can  obtain  the  relative  measure 
of  any  of  these  processes  with  respect  to  any  other  in  the  form  of 
nondlmensi onal  groups.  Spalding  [8]  identifies  105  such  groups  which  may  have 
an  apparent  relevance  in  a  general  combustion  scaling  problem.  It  is  fitting 
to  note  that  the  various  physical  and  chemical  processes  deemed  to  arise  in 
combustion  have  been  identified  intuitively  by  Spalding. 

F.  A.  Williams  [9],  in  contrast,  considers  the  mass  fire  problem  which  is 
somewhat  more  specific  in  many  respects  and  yet  more  complex  in  other  respects 
than  Spalding's  general  combustion  problem.  Rased  on  equations  of  change, 

Williams  identifies  a  total  of  29  nondimensional  groups,*  19  of  which  arise 
from  the  conservation  equations  themselves  and  10  from  the  boundary 
conditions.  (Of  these  last  10,  six  pertain  to  the  fuel  element  description 
and  four  to  the  ambient  atmosphere.) 

Since  these  many  identifiable  nondimensional  groups  far  surpass  the 
experimentally  plausible  degrees  of  freedom,  one  has  to  discern  the  scaling 
groups  which  need  not  be  preserved  invariant  between  the  model  and  the 

•Williams'  enumeration  in  fact  involves  far  more  than  29  groups,  the  actual 
number  being  determined  by  the  total  number  of  concerned  chemical  species. 
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prototype  experiments.  It  Is  in  this  discernment  that  one  needs  physical 
Insight  and  practical  experience.  It  Is  also  in  this  discernment  that  one 
recognizes  the  impossibility  of  complete  modeling  of  a  combustion  problem  and 
the  need  to  be  satisfied  with  'partial  modeling*. 

Even  partial  modeling  is  seldom  possible  of  the  entire  chronological 
behavior  of  a  fire  from  ignition  to  sustainment  to  spread  and  growth  to 
quasi-steady  burning  to  subsidence  and  extinguishment.  One  is  thus  forced  to 
concentrate  on  one  stage  of  the  fire  —  say,  steady  burning  —  or  one  aspect 
of  it  --  say,  the  burning  intensity  —  and  settle  to  model  it  partially.  It 
is  not  at  all  uncommon  to  discover  that  the  scaling  constraints  fulfilling  the 
requirements  of  similarity  between  the  model  and  prototype  with  respect,  to  one 
aspect  of  fire  behavior  may  conflict  with  restraints  to  obtain  similarity  of 
another  aspect.  For  example,  scaling  rules  to  reproduce  the  convective 
effects  are  generally  incompatible  with  those  to  preserve  the  thermal 
radiation  effects. 

Thus,  as  Spalding  concludes,  the  permissibility  of  a  particular  scaling 
technique  depends  for  the  most  part  on  the  stated  limited  objectives  of  the 
study  and  accuracy  sought. 

The  scaling  considerations  discussed  in  the  rest  of  this  report  are  based 
on  the  premise  of  partial  modeling.  As  such,  the  discussion  follows  the 
scaling  restraints  associated  with  various  chronological  aspects  of  the  fire 
in  the  respective  independent  sections. 

In  each  section,  the  most  significant,  existing  knowledge  is  identified 
and  this  knowledge  is  converted  into  a  scaling  exercise  to  identify  the 
similarity  groups  which  are  important  in  the  design  of  a  scale  model 
experiment . 
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SECTION  5 


IRRADIANCE  FROM  THE  FIREBALL  TO  THE  TARGET 

Thermal  radiation  from  the  air  burst  of  a  nuclear  weapon  comes  in  its 
characteristic  double-peak  (Glasstone/Dolan  [10])  to  heat  all  the  materials 
exposed  to  it.  It  is  well-known  (Martin  [11])  that  the  second  temporal  peak 
contributes  most  to  this  heating  process.  If  the  heated  materials  are 
combustible,  they  pyrolyze  and  spontaneously  ignite  within  the  range  where  the 
thermal  fluence  exceeds  about  12  Joules/cm2.  Delineation  of  this  range 
requires  an  examintion  of  three  aspects:  first,  the  characteristics  of  the 
fireball  as  a  source  of  thermal  radiation;  second,  transmission  of  the  emitted 
thermal  radiation  through  the  atmosphere  to  the  target  combustible  materials; 
and  third,  the  thermal  response  and  ignition  behavior  of  the  exposed 
materials.  In  this  chapter,  we  deal  with  the  first  two  of  these  three  aspects 
deferring  consideration  of  the  third  to  the  following  chapter. 

Let  W  be  the  yield  of  the  weapon  usually  expressed  in  units  of  kilo-ton 
(kt)  equivalent  of  TNT  explosive,*  The  fraction  of  W  which  is  emitted  as 
thermal  radiant  energy  Etotal  depends  not  only  on  W  Itself  but  also  on  the 
height  of  burst  H.  Denoting  this  fraction  by  f  (  =  Etotal/vO«  Glasstone  and 
Dolan  (Page  313)  present  f  =  f  (W,H)  in  tabular  form.  In  general,  f  is 
larger  if  H  and  W  are  larger. 

Of  importance  in  incindiary  effects  of  nuclear  explosions  is  the  manner 
in  which  thermal  energy  is  released  with  respect  to  time.  Experience 
Indicates  that  the  rate  of  thermal  energy  release  P  exhibits  a  maximum  at  a 
time  t  =  tmax  and  then  gradually  decays.  Roth  tmax  and  Pmax  depend  upon  the 
yield  U  and  the  height  of  burst  H.  Glasstone  and  Dolan  give  (Page  309)  this 
dependency  in  the  following  form. 

F  =  c(-e-  )n  Wm  (5.1) 

po 

where  F  can  be  either  Pmax  or  t^ax*  p  is  the  density  of  air  at  the  height  of 
burst  and  p0  is  that  at  the  sea  level.  The  constants  c,  n  and  m  depend  on 

*1  kt  TNT  equivalent  =  1012  calories  =  4.186  x  1012  Joules. 
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whether  Pmax  or  fmax  Is  under  consideration  and  also  on  whether  the  height  of 
burst  is  below  or  above  15,000  ft.  Table  1  gives  these  constants. 

The  instantaneous  thermal  power  output  P  (Watts)  of  the  fireball  of 
Instantaneous  temperature  T(K)  and  radius  R(m)  can  be  expected  to  be  given  by 

P  «  W4R2  (5,2> 

where  a  is  the  Stefan-Boltzmann  constant  (a  3  5.67  x  10‘8  W/m2*4).  As  T  and  R 
change  with  time  so  also  does  P.  Apparent  in  this  P(T,R)  relation  is  the 
initial  Increase  in  P  due  to  the  increase  in  both  T  and  R  and  the  later  decay 
due  to  the  decrease  in  T4  overwhelming  the  increase  in  R 2.  This  temporal 
distribution  is  given  by  a  single  curve  (see  Fig.  7.84  on  page  311  of 
Glasstone  and  Dolan)  in  the  coordinates  of  P/Pmax  versus  t/tmax  to  absorb  the 
dependence  of  history  on  W  and  H. 

The  integral  of  P  from  t  =  0  to  t  gives  the  thermal  energy  emitted  to  the 
time  t  of  consideration. 

/  P(t)dt  3  E(t)  (5,3) 

o 

E(«)  being  the  quantity  identified  earlier  by  the  symbol  Etotal  =  fw*  The 
amount  of  thermal  energy  released  to  time  t  expressed  as  a  percentage  of  the 
total  release,  i.e.,  E(t)/E(»),  is  also  given  in  the  figure  mentioned  above. 

A  more  explicit,  albeit  dimensional,  result  of  E(t)  for  air  bursts  of 
different  yields  (but  at  altitudes  of  bursts  below  100,000  ft)  is  presented  in 
Figure  7.87  on  page  314  of  Glasstone  and  Dolan.  Evident  from  this  figure  is 
the  dependence  of  tmax  on  W  which  translates  into  the  fact  that  explosions  of 
high  yield  W  result  in  longer  pulse  lengths. 

Thus  it  is  clear  that  once  the  yield  W  and  height  of  burst  H  are  given, 
one  can  read  off  the  thermal  partition  f  from  a  table  to  obtain  the  total 
thermal  energy  output  E(«0  as  fxW.  Simultaneously,  the  peak  power  Pmax  and 
the  time  to  its  occurrence  tmax  can  be  calculated  from  the  known  c,  n  and  m  of 
of  Table  I.  One  can  then  enter  into  Figure  7.84  of  Glasstone  and  Dolan  at  any 
time  t  to  read  off  P(t)/Pmax  as  well  as  E(t)/E(»).  Having  already  estimated 
Pmax  and  E(-),  the  instantaneous  power  P(t)  and  thermal  energy  release  E(t)  to 
time  t,  are  immediately  obtainable. 


If  the  fireball  Is  taken  as  a  point  source  emitter,  the  radiant  flux 
arriving  at  a  target  located  at  a  distance  0  away,*  Is  given  by 

I(t)  »  P(t)  x  t/4*D2  (5.4) 

where  x  is  the  transmissivity  of  the  intervening  air  mass  and  the  denominator 
stands  for  the  area  of  the  imaginary  spherical  surface.  This  relation  is 
exact  for  a  point  or  spherically  symmetric  source.  When  the  source  is  real, 

It  is  a  reasonable  approximation  for  large  values  of  0.  The  transmissivity  t 
of  the  atmosphere  is  less  than  unity  due  to  the  combined  effects  of  absorption 
and  scattering.  If  scattering  were  absent,  then  t  3  exp(-«D)  were  k  is  the 
absorption  coefficient;  this  corresponds  to  an  atmosphere  free  of  clouds  and 
moisture.  Figure  7.98  of  Glasstone  and  Dolan  shows  the  transmissivity  on  a 
clear  day  as  a  function  of  the  height  of  burst  and  the  distance  from  ground 
zero. 

The  units  of  I,  of  course,  are  W/m2.  From  the  known  P(t)  history  due  to 
a  brust  of  yield  W  at  a  given  altitude  H,  applying  the  -r(H,d)  from  Figure  7.98 
of  Glasstone  and  Dolan  and  knowing  D  =  /[H2  +  d2],  the  instantaneous 
irradiance  I  -  I(t;W,H,d)  is  thus  obtainable.** 

Once  I  is  mapped  out,  we  are  faced  with  the  response  of  target  materials 
exposed  to  it.  A  discussion  of  this  response  constitutes  the  content  of  the 
following  chapter.  From  laboratory  data  on  the  ignition  response  of  various 
materials,  Martin  developed  the  range  within  which  ignition  is  possible  from  a 
10  MT  air  bursts  at  different  altitudes.  This  result  is  presented  in  Figure  1 
on  page  3.17  of  Ref.  [11].  The  limitations  associated  with  this 
development  and  the  scaling  considerations  necessary  to  validate  it  with 
experiments  are  discussed  in  the  following  chapter. 


**  The  integral  I(t)dt  is  known  as  fluence,  F(t)  whose  units  are  J/m2. 
o 

This  quantity  physically  represents  the  total  energy  deposited,  from  time 
equal  to  0  to  t,  on  a  unit  area  of  the  target  surface.  The  fluence  is  to 
irradiance  as  E  is  to  oower.  Thus 
a  t  t 

F(t)  =  /  I (t)dt  =  /  P(t)  tMffDZdt  =  E(t)  t/4 tiO  2  . 
o  o 

With  E(t)  and  t  obtained  respectively  from  figures  7.84  and  7.98  of 
Glasstone  and  Dolan,  the  instantaneous  fluence  at  any  distance  d  from  a 
burst  of  a  given  yield  W  at  a  height  His  thus  readily  calculable.  The 
total  fluence  ?(*>)  is  simply  E(«)t/4ttD2  3  fWT/4^2. 
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It  Is  clear  from  the  preceding  discussion  that  one  can  make  a 
satisfactory  estimation  of  the  irradlance  due  to  a  single  burst  in  the  air. 
Existing  ability  to  estimate  the  irradlance  due  to  multiple  bursts  Is  not  as 
satisfactory.  The  alteration  of  the  atmospheric  attenuation  Is  perhaps  the 
most  difficult  issue  associated  with  this  multiburst  problem. 


SECTION  6 
IGNITION 


Spontaneous  Ignition  of  radiantly  heated  combustible  solids  has  been  an 
extensively  studied  subject  In  the  fields  of  both  propellant  combustion  and 
fire  safety.  Much  of  the  fire-related  ignition  study  has  been  of  experimental 
nature.  Certain  special  observations  have  to  be  enumerated  here. 

(a)  Shape  of  the  irradiance  pulse:  The  complicated  shape  of  thermal 
energy  release  from  an  air  blast  discussed  in  the  preceding  chapter  is 
difficult  to  reproduce  in  a  laboratory  ignition  test.  While  the  pulse  shape 
discussed  in  ENW  and  summarized  in  the  previous  chapter  of  this  report  is 
reasonable  at  low  altitudes,  the  pulse  may  be  quite  different  in  shape  If  the 
altitude  of  burst  is  high.  As  Martin  [11]  states,  it  is  clear  from  a  large 
body  of  experimental  work  that  the  response  of  materials  to  the  ENW  pulses  can 
be  correlated  by  techniques  which  apply  to  rectangular  (or  square-wave  or 
constant-flux)  pulses. 

(b)  First  ignited  materials:  The  extensive  experimental  studies 
reported  by  Martin  suggest  that  the  most  crucial  materials  of  concern  in 
ignition  by  the  fireball  are  thin  cellulosic  and  related  polymeric  kindling 
fuels.  These  are  generally  capable  of  igniting  within  a  few  seconds  of 
exposure  to  fluences  less  than  about  120  W/cm2. 

(c)  Main  controlling  process:  The  experiments  also  showed  that 
conduction  of  heat  Into  the  ignitible  solid  is  the  primary  process  controlling 
the  spontaneous  Ignition  due  to  radiant  exposure.  Three  distinct  kinds  of 
ignition  -transient  flaming,  persistent  flaming  and  thin-body  (glowing) 
ignitions  -are  identified  as  dependent  upon  the  thickness  a  and  conductivity 
Ks  of  the  specimen  and  the  irradiance  I  and  duration  t  of  the  exposure. 

Martin  [1]  empirically  delineates  the  conditions  under  which  the  different 
kinds  of  ignition  occur. 

(d)  The  prominence  of  transient  conduction  in  spontaneous  ignition  leads 
to  infer  that  ignition  is  incipient  when  the  exposed  surface  attains  a 
prescribed  ignition  temperature.  Rased  on  this  criterion,  the  present  author 
[2]  attempted  to  theoretical ly  deduce  Martin's  map.  The  nondimensional 
parameters  which  require  to  be  preserved  invariant  between  different 
situations  to  scale  the  ignition  process  are  evident  from  this  analysis.  At 
the  minimum,  they  are: 


Exposure  parameter 

U/KS(T*-T0) 

(6.1) 

Fluence  parameter 

I^t/KS(T*-T0), 

(6.2) 

Convective  heat  loss  parameter 

lu/Ks 

(6.3) 

Thermochemical  parameter 

<VTo,/(T*-To> 

(6.4) 

T*  is  the  prescribed  ignition  temperature,  T0  is  the  solid's  initial 
temperature  and  Tp  is  the  solid  temperature  at  which  its  pyrolysis  is  copious, 
h  is  the  convective  heat  loss  coefficient,  t  is  time,  a  is  thermal 
diffusivity.  The  subscript  s  stands  for  the  solid. 

Note  that  the  ratio  of  the  fluence  and  exposure  parameters  yields  the 
well-known  Fourier  modulus,  cist/*2,  i.e.,  the  nondimensional  time  as 
traditionally  defined  in  diffusion  problems.  The  exposure  parameter  is  a 
characteristic  Biot  number  which  indicates  the  magnitude  of  the  energy 
deposition  rate  I  relative  to  the  characteristic  energy  conduction  rate  into 
the  solid  Ks(T*-T0)/z.  A  physically  thin  solid  of  high  conductivity  and 
characteristic  temperature  difference  (T*-T0),  when  exposed  to  a  low  intensity 
of  irradiation,  results  in  a  small  exposure  parameter.  The  solid  then 
entertains  negligible  internal  temperature  gradients  and  is  called  'thermally 
thin'.  A  physically  thick  solid  of  low  thermal  conductivity  when  heated 
intensely,  results  ir.  a  large  value  of  exposure  parameter.  Internal 
temperature  gradients  predominate  in  such  a  'thermally  thick'  solid. 

Upon  incorporating  the  definition  of  thermal  diffusivity  eKs/psCpS, 
(where  the  denominator  is  the  volumetric  heat  capacity  of  the  solid),  one  can 
see  the  fluence  parameter  take  the  form  It/psCps(T*»T0)  t.  In  this  form,  it  is 
evident  that  the  fluence  parameter  is  the  ratio  of  total  energy  supplied  to 
the  solid  to  the  energy  involved  in  raising  the  solid  temperature  from  T0  to 
T*. 

The  convective  loss  parameter  apparently  becomes  important  only  when  the 
solid  is  thermally  thin.  It  obviously  depends  upon  not  only  the  orientation 
of  the  sample  with  respect  to  gravity  but  also  on  the  characteri stic 
dimension.  The  gas  phase  thermophysical  properties  also  intrude  into  the  heat 
transfer  loss  coefficient  h.  It  is  interesting  to  note  that  the  ratio  of  the 
loss  parameter  to  the  exposure  parameter  yields  h(T*-T0)/I  which  represents 
the  rate  of  energy  loss  by  convection 
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(characteristically  when  T=T*)  expressed  as  a  fraction  of  the  irradiant  energy 
supply  rate. 

If  only  a  restricted  class  of  materials  is  of  concern,  then  the 
thermochemical  parameter  can  be  taken  as  fixed.  There  are,  however, 
approximations  pertaining  to  the  pyrolysis  kinetics  (via  Tp)  and  to  the 
kinetics  of  oxidation  of  the  pyrolyzate  efflux  (via  T*)  embedded  in  this 
parameter. 

(e)  The  experience  of  difficulty  in  igniting  fuels  which  are  thick  is 
evidently  explained  by  the  excessive  conductive  drain  into  the  solid.  That 
moist  fuels  are  harder  to  ignite  is  also  perhaps  explainable  through  the 
influence  of  moisture  on  T*  and/or  Tp. 

(f)  Even  more  important  is  the  implication  of  points  (c)  and  (d)  that  the 
fluence  required  to  produce  ignition  is  higher  if  the  exposure  flux  is  lower. 
In  fact,  there  exists  enough  convincing  experimental  evidence  to  indicate  that 
the  ignition  temperature  T*  is  not  altogether  a  constant  for  a  given  material 
as  embodied  the  delineations  of  Refs.  [1  and  2];  and  that  it  is  higher  at  a 
lower  irradiance,  such  a  dependency  being  more  pronounced  at  smaller 
irradiances. 

(g)  A  serious  consequence  of  this  T*(I)  dependency  is  the  existence  of  a 
threshold  irradiance  level  below  which  ignition  of  the  solid  v;ill  not  occur 
even  with  an  indefinitely  long  duration  of  exposure. 

(h)  The  issues  of  Ignition  thus  are  two-fold.  First,  does  a  given 
exposure  irradiance  produce  ignition  of  a  given  target?  Second,  if  the  answer 
to  the  first  question  is  in  the  affirmative,  how  long  an  exposure  duration  is 
required  to  produce  ignition?  While  the  first  question  is  concerned  with  the 
ignition  threshold,  the  second  is  with  the  ignition  time. 

A  recent  work  by  Gandhi  [3]  addresses  these  two  issues.  One  of  the 
accomplishments  if  this  work  is  the  description  of  the  energy  balance 
associated  with  heating  the  solid  radiatively  to  ignition.  As  the  solid  gets 
heated  by  transient  conduction,  its  hot  surface  loses  energy  not  only  by 
thermal  reradiation  to  the  surroundings  but  also  by  natural  convection  of  the 
adjacent  gas.  Among  the  primary  descriptors  of  the  problem  are:  (a)  the 
orientation  of  the  element  which  influences  the  inception  and  effectiveness  of 
the  free  convective  loss  as  well  as  the  angle  of  impingement  of  the  radiant 
beam;  (b)  the  exposure  intensity  I  nondimensional ized  (with  the  element 
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conductivity  Ks  and  ambient  temperature  TM)  to  the  form  Iz/KsT*;  (c)  the  ; 

length  of  the  element  H  nondimensional 1  zed  (with  the  gravity  constant  g  and  ! 

gas-phase  kinematic  viscosity  vg)  to  form  a  Galileo  number  gHV\$3;  (d)  the  ; 

surface  emissivity  es*  normalized  (with  the  Stefan-Boltzmann  constant  a,  ; 

ambient  temperature  T„,  element  thickness  i  and  solid  conductivity  Ks)  to  form  i 

ae^TH/Ks*  T i me  t  is  normalized  to  be  a  Fourier  modulus  t/z2  where  is  ! 

the  virgin  solid  thermal  diffusivlty.  There  are  six  additional  parameters  of 
secondary  importance  contributed  by  the  kinetics  and  energetics  of  pyrolysis  ; 

and  oxidation  reactions.  The  ratios  of  various  thermophysical  properties  of  l 

the  solid  and  gas  phases  constitute  tertiary  parameters.  While  a  number  of  | 

important  predictions  are  possible  from  Gandhi's  theory,  it  is  sufficient  to 
note  here  that  the  most  crucial  v  "ing  parameters  to  study  ignition  are  the 
following. 

Orientation  ---  ! 


Exposure  parameter 

u/u 

s  =0 

(6.5) 

Height  parameter 

gH  v 

(6.6) 

Reradiant  loss  parameter 

<KT^Ks 

(6.7) 

Thickness  parameter 

Z/H 

(6.8) 

Time 

c^t/z2 

(6.9) 

(i)  The  orientation  of  the  target  element  is  known  to  influence  the 
possibility  of  ignition.  Martin,  for  example,  reports  on  experiments  in  which 
a  given  irradiance  may  be  capable  of  producing  ignition  of  a  larger  element  in 
vertical  orientation  but  may  merely  char  without  ignition  when  the  element  is 
horizontal  or  short. 

(j)  There  exist  a  number  of  field  tests  in  which  the  vulnerability  of 
combustibles  exposed  to  fireball  radiation  is  monitored  by  employing  beds  of 
shredded  paper,  shaved  wood  and  excelsior  as  targets.  In  view  of  the  preceding 
discussion,  the  factors  of  concern  in  this  practice  include:  (i)  the  bed 
orientation;  (ii)  the  packing  density  of  the  bed;  (iii)  the  bed  thickness,  if 
packing  is  dense;  (iv)  the  fineness  of  shredding;  and  (v)  the  paper  thickness. 
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It  Is  obvious  that  two  conductivities  of  the  condensed  phase  become  relevant; 
first,  the  effective  conductivity  of  the  porous  bed  and  second,  the 
conductivity  of  the  shred  element  material.  If  packing  is  rather  dense,  it  is 
the  overall  bed  characteristics  which  govern  the  ignition  process.  If  the 
packing  is  loose,  on  the  other  hand,  it  is  the  characteristics  of  the  shred 
element  which  are  critical.  Design  of  these  target  beds  to  maximize  the 
obtainable  information  is  by  Itself  an  important  and  interesting  research 
project. 
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SECTION  7 


IMPACT  OF  THE  BLAST  WAVE  ON  THE  NASCENT  FIRES 

Following  the  radiant  ignition,  the  scaling  of  which  has  been  discussed 
in  the  proceeding  chapter,  the  fires  would  grow.  (Fire  propagation  will  be 
discussed  in  the  following  Section.  The  blast  wave,  along  with  its 
associated  transient  pressure  disturbance  and  winds,  would  soon  arrive  to 
Impinge  upon  the  flames  which  are  but  only  in  their  earliest  stages  of  growth. 
Whether  or  not  the  blast  impingement  would  extinguish  the  flames  is  a  question 
of  concern. 

The  properties  of  concern  in  the  blast/fire  interaction  process  may  be 
conveniently  placed  into  three  categories  respectively  pertaining  to  the  fuel, 
fire,  and  the  blast.  Among  the  fuel  propeties  are:  the  pyrolysis  kinetics  of 
wood  fuel  or  the  vaporization  thermodynamics  of  liquid  fuels,  thermodynamics 
of  mixing,  thermochemi stry  of  glowing  and  flaming  combustion,  oxidation 
kinetics  of  glowing  and  flaming  reactions,  and  transport  properties.  The  fire 
properties  mainly  concern  the  geometry  (shape,  size,  orientation)  of  the  fuel 
bed,  the  nature  of  flow  (forced  or  free  convection)  and  the  externally  imposed 
radiant  flux.  The  blast  properties  include  the  peak  overpressure,  time  of 
positive  phase  and  continued  radiant  pulse.  Naturally,  the  consequence  of 
fuel  beds  redistributed  into  debris  piles  in  a  multiple  blast  situation  is  to 
be  considered  after  some  degree  of  understanding  is  obtained  about  the 
well-defined  single  fire  response  to  a  well-defined  single  blast  impingement. 

The  most  abundant  and  relevant  experiments  to  study  the  blast/fire 
interaction  process  are  those  reported  by  Martin,  Backovsky  and  their 
colleagues  at  SRI  International  [12-14],  The  SRI  experimental  facility 
attempts  to  capture  the  essential  variables;  the  air-driven  shocktube  is 
capable  of  producing  shockwaves  typical  of  those  produced  by  ki 1 o-to-megaton 
nuclear  explosions  in  air.  The  peak  over  pressures  (up  to  25  psi)  and 
positive  phase  durations  (between  0.10  to  over  3.5  sec)  are  prescribable  by 
the  operator  in  this  experimental  set  up.  Liquid  pool  fires  are  tested  with  a 
variety  of  fuels:  n-hexane,  methanol,  kerosene,  n-pentane,  and  acetone  with 
the  fuel  bed  length  ranging  between  12  and  36  inches.  The  mean  overpressure 
is  varied  between  0.91  and  7.5  psi.  Time  of  positive  phase  duration  was 
ranged  between  0.068  and  3.8  sec.  In  some  tests,  the  preburn  time  before  the 
blast  arrival  is  noted.  Barriers  located  upstream  of  the  pool  are  also 
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studied  in  some  tests.  The  observations  included  noting  whether  or  not  a 
given  fire  is  extinguished  by  a  given  blast. 

Fires  of  typical  wood  cribs  made  of  (mostly)  dry  western  hemlock  and 
(some  redwood  are  subjected  to  overpressures  in  the  range  0.91  to  9.88  psi 
with  postive  phases  between  0.078-0.117  sec.  The  cribs  are  made  of  sticks 
either  3/8"  or  3/4"  thick.  A  tray  of  propyl  alcohol  under  the  crib  Is  used  to 
prompt  the  crib  ignition  following  which  the  crib  is  allowed  to  burn  for 
anywhere  from  about  60  to  180  sec  before  the  blast  wave  is  imposed  on  the 
fire. 

There  are  also  reprted  eleven  tests  in  which  trays  of  shredded  blotter 
or  filter  paper  act  as  fuel  beds. 

In  an  attempt  to  identify  the  mechanisms  of  interaction  between  a  fire 
and  a  blast,  Kanury  [5]  developed  a  preliminary  scaling  analysis  in  which  the 
physics  and  chemistry  of  combustion,  as  perturbed  by  the  blast  effects,  come 
to  Interplay  and  decide  the  fate  of  the  fire.  The  blast  brings  with  it  a  rise 
in  pressure  and  temperature  as  well  as  winds.  While  the  pressure  and 
temperature  exert  an  influence  on  the  oxidation  reaction  rate  in  a  more  or 
less  direct  way,  the  air  flow  exerts  a  more  complicated  influence  on  both  the 
chemistry  and  physics.  F.nhanced  convection  dilutes  the  reactant  mixture 
thermally  as  well  as  compositionally.  Reduction  in  temperature  is  the 
consequence  of  thermal  dilution.  Increase  in  the  concentration  of  oxygen  (of 
air)  is  the  result  of  the  compositional  dilution.  These  dilutions  directly 
reflect  in  conditions  conducive  to  the  extinguishment  of  a  flame  (namely,  the 
gas-phase  oxidation  reaction).  However,  when  the  fuel  bed  is  constituted  of 
such  materials  as  wood  whose  char  may  undergo  glowing  combustion,  the 
augmented  oxygen  supply  tends  to  augment  the  combustion  intensity.  (This 
phenomenon  is  quite  familiar  to  all  every-day  combustion  experimental i sts  who 
play  with  the  barbeque  charcoal  grills.)  The  outcome  of  these 
physico-chemical  interactions  of  the  blast  with  the  fire  is:  (a)  the  blast, 
if  sufficiently  strong,  would  blow  away  the  flames;  and  (b)  it  would 
exasperate  the  glowing  process.  In  fact,  if  the  augmentation  of  glowing  is 
sufficiently  large,  the  glowing  objects  would  support  reflashing  of  a  flame 
upon  subsidence  of  the  blast  effects.  Based  on  a  scrutiny  of  this  interplay 
between  physical  transport  processes  and  chemical  kinetic  reaction  processes, 
Kanury  [6]  identifies  a  Damkohler  number  P  to  represent  the  inverse  of  the 
blast  strength  and  a  fire  strength  parameter  q  which  defines  the 
energeticity  of  the  fire.  The  definitions  of  these  two  scaling  parameters  are 
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as  following: 


Inverse  blast  strength  P  h  k0jj/u  (7.1) 

Fire  strength  q  =  RhcYp-( /ECpg  (7.2) 

k0  Is  the  preexponential  factor  of  an  assumed  first-order,  overall  oxidation 
reaction,  i  is  the  fuel  bed  characteri stic  dimension  and  u  is  the  maximum 

velocity  of  wind  induced  by  the  blast  wave.  R  is  the  universal  gas  constant, 
hc  is  enthalpy  of  combustion,  is  the  fuel  mass  fraction  at  the  surface,  E 

is  the  activation  energy  of  the  oxidation  reaction  and  Cpg  is  the  gas-phase 

specific  heat.  As  given,  P  is  a  Damkohler  numer,  the  ratio  of  chemical  to 

physical  speed,  u  will  depend  upon  the  pressure  as  perturbed  by  the  blast. 

On  a  map  of  P(q),  regimes  in  which  extinction  is  expected  are 
distinguished  from  those  in  which  the  fire  will  sustain  the  impact  of  the 
blast.  Where  charring  fuels  are  involved,  the  role  of  glowing  combustion  in 
the  reflashing  of  an  extinguished  flame  is  examined.  The  following 
conclusions,  which  are  reasonably  well-founded  by  the  SRI  experimental 
observations,  are  reached. 

(a)  Flames,  on  both  liquid  and  solid  fuels,  appear  to  be  blown  out  by  the 
same  mechanism. 

(b)  The  P  required  to  achieve  extinction  is  lower  for  flames  of  higher  q. 

(c)  The  influence  of  an  upstream  barricade  is  shown  to  render  the  liquid 
fuel  flames  more  blast-resistant  so  that  the  nonextinction  region  of  the  P(q) 
map  is  increased.  No  distinct  boundaries,  however,  can  now  be  drawn  for  this 
barricade  effect.  Experiments  are  required  in  which  taller  barricades  are 
systematically  tested.  Theory  describing  the  mechanisms  of  flame 
stabilization  by  reci rculation  has  to  be  developed  to  compare  with  these 
experiments . 

(d)  The  enthalpy  of  combustion  of  the  fuel,  fuel  bed  length,  barricades 
and  shock  wave  strength  are  the  most  important  variables  in  the  blast/class  B 
fire  interaction  process.  The  fuel  vaporization  thermodynamics,  including 
enthalpy  of  vaporization,  the  nature  of  wick  holding  the  fuel  in  the  pool,  and 
combustion  kinetic  parameters  appear  to  have  little  or  no  effect  on  the 
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response  of  class  8  fires  to  blasts.  This  is  a  significant  point  since  such 
properties  as  those  pertaining  to  global  oxidation  kinetics  can  not  be 
expected  to  be  known  in  reliable  and  quantitative  detail  ever. 

(e)  The  hypothesis  of  blast/fire  interaction  postulated  and  tested  with 
pool  fire  data  is  capable  of  describing  fire  extinguishment  viewed  in 
different  ways,  for  instance,  if  extinguishment  is  defined  as  annihilation  of 
the  flame  it  Includes  both  the  degeneration  of  a  stationary  flame  and 

' nonreappearance'  of  a  displaced  flame.  The  flame  displacement  mechanism 
advocated  by  the  SRI  group  thus  is  completely  consistent  with  the  hypothesis. 

(f)  Having  correlated  the  12  to  36  inch  long  liquid  pool  fire  in  Ref. 

[5],  the  next  task  in  studying  class  B  fires  Is  to  develop  ways  of  scaling  the 

fires  up  to  several  tens  or  hundreds  of  feet  in  size.  Development  of  these 

scaling  rules  requires  consideration  of  such  additional  factors  as  :  thermal 

radiation  feed-back  from  the  larger  flames  to  the  fuel  surface,  effect  of  the 
size  of  the  flames  on  the  uniformity  with  which  the  blast  would  perturb  them, 
accentuation  of  the  recirculation  phenomena,  and  others.  Based  on  the 
definition  of  P,  the  critical  fuel  bed  length  is  nearly  proportional  to  the 
peak  over  pressure  and  the  particle  velocity.  This  conclusion  is  similar  to 
that  reached  earlier  by  Backovsky  et  al  in  Ref.  [131. 

(g)  Turning  to  class  A  fires  as  represented  by  cribs  of  wood  sticks,  It 
is  noted  that  the  mechanisms  of  fl ame  extinction  appears  to  be  the  same  as 
that  of  class  B  fires.  The  extinguished  state,  however,  involves  charcoal 
embers  capable  of  undergoing  a  metastable  glowing  combustion  process.  This 
glowing  combustion  process  supplies  sufficient  energy  to  the  solid  to  continue 
production  of  pyrolyzates.  As  the  blast  effects  subside,  the  cooling  (and  yet 
pyrolyzing)  embers  may  pass  through  a  thermochemical  state  corresponding  to 
reflashing  of  flame  in  the  mixture  of  pyrolyzate  and  air. 

Based  on  the  knowledge  of  the  composition  of  pyrolysis  products,  and  of 
the  relative  independence  of  the  extinction  boundary  of  class  B  fires  on  the 
specific  nature  of  the  fuel,  there  seems  to  be  no  compelling  reason  to  believe 
that  the  flame  extinction  boundary  for  crib  fires  will  be  different  from  that 
for  class  B  fires.  The  correlation  of  P  versus  q  confirms  this  idea  although 
this  correlation  is  not  as  clear-cut  as  that  of  liquid  fuel  flames.  The 
reasons  for  this  can  be  any  combination  of  such  complications  as  the 
following,  (i)  Although  it  appears  logical  to  expect  the  preburn  time  t.o  be 
an  indirect  measure  of  the  fire  strength  q,  the  exact  connection  is  not  clear 


22 


at  present.  Further  study  Is  required  of  this  Issue  to  establish  the  role  of 
transient  conduction  (with  pyrolysis)  in  rendering  persistence  and  strength  to 
a  crib  fire.  (11)  In  a  crib,  opportunity  exists  for  the  flames  to  find 
intense  recirculation  zones  behind  individual  sticks.  These  reci rctil ation 
zones  are  more  effective  if  the  sticks  are  thicker.  Because  recirculation 
tends  to  stabilize  the  flames  and  thus  to  shift  the  extinction  boundary,  it  is 
reasonable  to  expect  nonextinctions  where  extinction  would  be  expected  in  the 
absence  of  reci rculation.  The  future  study  recommended  above,  of  flame 
stabilization  by  reel rcul ati on ,  will  aid  in  understanding  the  effect  of  stick 
thickness  on  extinction.  The  high  surface  temperature  of  the  wood  sticks  in  a 
crib  will  be  an  important  additional  parameter  to  be  considered  in  this 
connection,  (ill)  The  surfaces  at  different  locations  in  the  crib  would  char 
nonuni formly  In  all  crib  fires.  The  distribution  of  glowing  patches  which  may 
expediently  pilot  the  reignition  of  blown-out  flames  is  thus  a  serious 
ambiguity  of  all  crib  fires.  Resolution  of  this  ambiguity  requires  further 
work,  (iv)  Compared  to  liquid  pools,  the  crib  fuel  beds  are  volumetric  by  two 
different  counts,  first  by  the  crib  structure  itself  and  second  by  the 
conduction-controlled  interior  of  the  individual  sticks.  While  the  flames 
over  pools  are  all-important  in  determining  the  pool  burning  rate,  the  flame 
over  the  crib  (i.e.,  the  'top-flame')  is  of  little  or  no  significance  compared 
to  the  flames  within  the  crib,  (v)  Thermal  radiation  from  the  flames  to  the 
fuel  surface  is  more  or  less  unimportant  in  the  class  B  fires  in  the  size 
range  tested.  In  contrast,  thermal  radiation  plays  a  dominant  role  in  the 
life,  strength  and  death  of  a  crib  fire.  Additional  study  is  required  to  gain 
a  better  understanding  of  this  role. 

The  data  correlations  of  Ref.  [5]  also  point  out  the  possibility  that  the 
basic  burning  process  of  shredded  paper  fuel  beds  may  be  altogether  different 
from  that  of  cribs.  It  may  be  closer  to  that  of  liquid  fuels,  but  this  is  yet 
to  be  examined  by  additional  study. 

(h)  The  flame  reflashing  mechanism  postulated  in  Ref.  [5]  can  not  at 
present  be  tested.  It  is  surprising,  however,  to  note  how  few  experiments  in 
the  SRI  collection  in  fact  show  a  distinct  reflashing.  In  as  much  of  charcoal 
glowing  is  expected  to  play  an  Important  role  in  this  reflashing  process,  it 
appears  important  that  a  detailed  study  be  undertaken  of  the  effect,  of  blast 
on  glowing  wood  surfaces.  The  nature  of  glowing  as  it  is  disturbed  by  the 
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last  variables,  the  implications  in  continued  pyrolysis,  and  possible 
transition  of  glowing  to  flaming,  are  not  only  important  in  uderstanding  the 
reflashing  process  but  also  exciting  from  the  view-point  of  science. 

(1)  Turning  to  scaling,  not  withstanding  the  complicated  nature  of  wood 
crib  fire  burning  discussed  above,  crib  fires  are  probably  easier  to  be  scaled 
up  since  their  burning  Is,  more  or  less,  locally  controlled.  Oxygen  depletion 
in  the  depths  of  the  crib  may,  however,  become  a  crucial  factor  in  such 
scaling  endeavors.  Thermal  radiation  effects,  of  course,  are  critically 
required.  One  important  conclusion  of  the  present  study  is  that  wood  crib 
burning  with  and  without  blast  interaction  requires  much  further  study  before 
results  useful  in  efforts  of  mitigation  can  be  obtained. 

(j)  Shredded  paper  tray  fires  are  in  a  class  all  by  themselves  in  their 
behavior  unexplained  at  present. 

The  conclusion  one  can  reach  from  this  status  report  is  that  while  some 
preliminary  aspects  of  blast/fire  interaction  have  been  studied  to  date,  much 
remains  to  be  done.  The  blast  strength  parameter  P,  fire  strength  parameter 
q,  and  barricade  height  (ratioed  with  the  fuel  bed  length)  appear  to  be  the 
predominant  scale  variables.  Shredded  paper  tray  fires  are  altogether  poorly 
understood  as  to  their  response  to  an  imposed  blast  wave.  Lastly,  little  is 
known  of  the  response  of  vertical  fuel  bed  fires  and  of  the  role  played  by 
thermal  radl ation--both  from  the  thermal  pulse  residual  and  from  the  fire 
itself.  Multiple  blast  effects  and  multiple  fire  effects  are  also  yet  to  be 
studied. 


SECTION  8 


FIRE  PROPAGATION  AND  GROWTH 

Born  of  the  Impingement  of  the  thermal  pulse  on  combustibles  and 
perturbed  in  its  infancy  by  the  arrival  of  the  blast  wave  with  its  attendant 
rise  in  pressure,  temperature  and  wind,  the  fire  finds  for  its  survival  a 
scenario  of  disturbed  objects  amidst  piles  of  debris.  The  same  blast  which 
caused  mechanical  disruption  of  bodies  to  generate  the  debris,  also  brings 
about  the  secondary  fires.  Not  withstanding  the  fact  that  combustibles  and 
noncombustibles  are  randomly  mixed  to  form  the  debris  piles,  these  piles  are 
expected  to  offer  an  effective  base  upon  which  both  the  secondary  fires  and 
the  blast-survived  primary  fires  flourish.  A  combination  of  both  smoldering 
and  flaming  fires  is  foreseeable,  since  the  fires  are  expected  to  continue 
burning  for  days  (or  even  weeks)  after  the  blast,  they  pose  a  continued  (and 
growing)  hazard  which  can  be  either  further  aggrevated  or  brought  under 
control . 

Fire  growth  is  a  topic  of  extensive  scientific  study.  The  mechanisms  of 
flame  spread  on  contiguous  fuel  beds  in  small  scale  have  been  thoroughly 
reviewed  by  F.  A.  Williams  [15].  The  relatively  better  definable  and 
controllable  configurations  of  downward  flame  spread  (on  vertical  combustible 
surfaces)  and  horizontal  spread  (on  vertical  or  horizontal  surfaces)  have  been 
studied  most  rigorously.  The  primary  conclusion  arrived  in  these  studies 
relates  to  the  mechanism  of  heat  transfer  from  the  'burning'  zone  to  the  fuel 
yet  to  be  Involved.  Conduction  through  the  gas  phase,  througri  the  solid 
phase,  convection  in  the  gas  phase,  convection  in  the  condensed  phase  (in 
liquids)  and  thermal  radiation  are  among  the  identified  hert  transfer 
mechanisms.  While  these  studies  are  useful  In  understanding  tne  fundamental 
processes  of  fire  growth,  they  are  less  relevant  in  practice  than  situations 
involving:  (a)  spread  of  fire  by  a  smolder  front;  (b)  horizontal  and  downward 
flame  spread  in  opposing  wind;  (c)  upward  fire  spread  on  vertical  surfaces; 

(d)  fire  spread  in  a  forest;  (e)  radiant  heating  and  jumping  fire  spread  in 
descrete  fuel  elements;  and  (f)  growth  in  compartments. 

8.1  Spread  of  a  Fire  by  a  Smolder  Front:  A  nonflaming  exothermic 

oxidation  reaction  wave  propagating  within  a  porous  combustible  solid  medium 
constitutes  smolder.  As  energy  is  released  at  the  smolder  front,  part  of  it 
is  lost,  predominantly  by  conduction,  through  the  solid  to  the  surroundings 
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while  the  rest  is  employed  to  preheat  the  solid  to  the  state  of  smoldering. 
Ohlemiller  and  Rogers  [16]  propose  that  if  the  balance  between  energy 
generation  and  loss  were  to  be  the  controlling  factor  in  determining  a 
sustained  smolder  spread,  then  the  energetics,  transport  properties  and 
kinetics  of  smolder  reactions  have  to  be  such  that  a  nondimenslonal  parameter 
5,  defined  below,  has  to  exceed  a  certain  minimum  value. 

h  ZE*2exp(-E/RT  ) 

.  _  0  P  to  1 


The  larger  the  heat  loss  (expressed  as  a  Biot  number  defined  below),  the 
larger  the  minimum  value  of  6. 

m  •  jji  .  <8'2: 

hc  Is  heat  of  combustion  of  the  smolder  reaction  whose  preexponential  factor 
is  Z  and  activation  energy  is  E.  R  is  gas  constant.  Tp  is  smolder 
temperature,  i  is  half-thickness  of  the  solid  whose  conductivity  is  K.  The 
surface  of  the  solid  loses  heat  to  the  surrounding  by  convection,  h  being  the 
corresponding  heat  transfer  coefficient.  A  high  ambient  temperature  T,,, 
nondimenslonal ized  as  9„  a  E(T„  -  Tp)/RTp2  causes  the  critical  s  to  he 
smaller.  In  contrast.  Palmer  [17]  and  other  workers  invoke  the  rate  of  oxygen 
transfer  through  the  porous  layer  as  a  governing  factor  of  the  smoldering 
rate.  Two  situations  can  be  envisioned.  In  the  frst,  the  smolder  wave 
propagates  into  the  material  from  near  the  surface  exposed  to  the  environment. 
In  this  case,  oxygen  diffusion  has  to  occur  through  the  residual  'char'  to  the 
propagating  front.  Since  the  char  thickness  offers  a  gradually  increasing 
resistance  to  oxygen  diffusion,  the  smolder  speed  decreases  as  the  Inverse 

square-root  of  time  and  hence  the  smoldered  thickness  varies  as  /t.  In  the 
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second,  the  smolder  propagates  from  within  the  fuel  bed  towards  the  exposed 
surface,  thus  encountering  progress! vely  smaller  diffusion  resistance  to  the 

oxygen  supply.  Although  a  /t  dependency  arises  again  (as  characteristic  of 
all  diffusion  problems),  the  spread  speed  increases  with  time  in  this 
situation.  The  smolder  spread  history  thus  appears  to  be  scaleable  in  terms 
of  the  dependent  parameter  V  and  independent  parameter  D  defined  below. 
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where  ps  and  pc  respectively  are  the  initial  and  final  density  of  the  solid; 
pg  is  the  ambient  gas  (i.e.,  air)  density;  ts  is  the  physical  characteristic 
thickness  of  the  smoldering  solid  of  specific  heat  Cps  and  conductivity  Ks;  D 
is  diffusivity  of  O2  through  the  char  or  the  original  solid  depending  upon  the 
spread  direction,  Y0o,  is  the  ambient  oxygen  mass  fraction  (normally  0.23),  and 
f  is  the  stoichiometric  mass  ratio  of  fuel  to  oxygen,  ti  is  spread  speed  and  t 
is  time. 

The  state-of-the-art  of  smolder  spread  is  still  in  its  infancy.  Detailed 
scaling  rules  are  yet  to  evolve,  especially  when  thermal  radiative  heat 
transfer  within  the  porous  fuel  bed  assumes  the  controlling  responsibility  for 
the  spread  process. 

8.2  Horizontal  and  Downward  Flame  Spread  in  Opposing  Mind:  As  a  prelude 

to  the  problems  of  more  direct  concern  in  the  context  of  this  report,  we  shall 
review  here  but  briefly  the  salient  current  understanding  of  the  wind-opposed 
fire  spread.  The  literature  abounds  with  reports  on  the  mechanism  of  this 
fire  spread  process. 

de  Ris  [18]  presented  a  landmark  analysis  of  fire  spread  on  a  horizontal 
fuel  bed  (thin  and  thick)  in  opposed  oxidant  flow  to  bring  the  first  orderly 


quantitative  interpretation  of  all  the  then  existing  flame  spread  information. 
The  primary  accomplishment  of  this  work  is  the  derivation  of 
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for  thin  fuels  of  thickness  zs.  The  evaporating  surface  temperature  Tvap  is 
representative  of  the  evaporation  thermodynamics  or  pyrolytic  decomposition. 
Tg  is  the  temperature  of  the  spreading  flame.  T*  is  temperature  of  the 
ambient  gas  whose  conductivity  is  Kg.  For  thick  fuel  beds,  the  following 
prediction  is  made. 


»sc.sKsv“  .  j  Tf  -  Tvap  .  281  *  2V»  j  (8.6: 

pgcpgKgug  Tvap  ~  T«  °gcpgug^Tvap  TJ  ’ 

where  subscript  g  stands  for  the  gas  phase,  Ug  is  the  velocity  of  opposing 
wind  flow,  Ri  and  R2  are  upstream  and  downstream  radiative  heat  fluxes,  and 
F(z)  is  a  given  function  of 

z  = _ _  (8.7! 

pC  u  L 

Mg  pg  g  1 

where  ai  is  the  characteristic  length  of  the  forward  radiation. 

These  remarkable  relations  vividly  point  out  the  manner  in  which  various 
properties,  conditions  and  variables  combine  to  form  the  scale  parameters  for 
the  phenomenon  of  laminar  horizontal  diffusion  flame  spread. 

F.  A.  Williams  [15]  pulled  together  the  various  mechanisms  of  fire  spread 
over  condensed-phase  combustibles  into  a  coherent  framework  based  on  the 
following  nondimensi onal  rendition  of  the  classical  concept  of  Mallard  and 
Le  Chatelier. 


psuahs 
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Here,  ps  is  the  fuel  density,  u  is  fire  spread  velocity,  q"  is  the  rate  of 
heat  transfer  across  the  surface  (or  line)  of  fire  inception  and  &hs  is  the 
rise  in  fuel  enthalpy  required  to  bring  the  supply  fuel  to  the  ignition  state. 


The  heat  flux  q"  depends  upon  the  specific  mechanism  of  heat  transfer 
which  drives  the  flame  spread  process.  If  it  is  the  solid-phase  conduction, 
then  q"  «  Ks(TVap  -  T.)/As  where  As  is  the  conduction  characteristic  length. 

It  is  thus  evident  that  psu  Ahs As/Ks(Tvap  -  Tj  -1  becomes  the  scaling 
prescription  in  this  situation.  When  Ahs  can  be  written  as  a  product  of 
specific  heat  Cps  and  temperature  rise,  this  prescription  further  reduces  to  a 
Peclet  number  criterion,  viz:  uAs/as  •  1  where  <*  3  Ks/pscps  is  thermal 
diffusivity.  The  conduction  length  As  Itself  may  either  be  a  geometric 
parameter  such  as  the  fuel  element  thickness  or  be  a  transient  parameter, 

i.e..  As  "  /«st,  where  t  is  time. 

If  conduction  in  the  gas  phase  were  to  be  operative,  the  heat  flux  is 
given  by  q"  »  h(T^.  -  Tyap)  so  that  ps  u  Ahs/h(T^.  -  Tyap)  «•  1  constitutes  the 
scaling  restraint.  Denoting  the  gas  phase  conduction  length  by  Ag>  it  is 
noted  that  h  *  Kg/Ag  where  Kg  is  the  gas  phase  conductivity  so  that  one  can 
arrive  at  the  alternate  scaling  law  ps  u  Ah$  Ag/Kg(Tf  -  Tvap)  «  1.  For  a 
thermally  thin  fuel  element,  Ag  may  be  expected  to  be  of  the  same  order  of 
magnitude  as  As  so  that  the  present  results  are  consistent  with  the  result  of 
de  Ris. 

Let  as  be  the  depth  to  which  a  thick  solid  (of  density  ps,  specific  heat 

Cps  and  conductivity  Ks)  required  to  be  heated  (from  TM  to  Tvap)  in  order  to 

receive  the  leading  edge  of  the  flame  spreading  at  a  speed  of  u.  Then  the 

energy  demand,  per  unit  width  of  the  solid,  is  given  by  ^  u  Cps  As(Tvap  “  Tj 

Watts  per  meter.  The  energy  supply  rate,  per  unit  width,  from  the  flame  to 

the  solid  at  the  leading  edge  is  d"g£g  (W/m)  where  the  heat  flux  q"g  will  be 

further  discussed  later;  and  the  characteristic  gas  phase  length  is  Ag.  The 

ratio  of  demand  to  supply  constitutes  a  dimensionless  ratio  of  importance  in 

the  spread  process.  This  ratio  is  p^  C  u  (T  -  T  )/q"  a_.  The 

s  ps  s  vap  *  9  9 


perturbed  solid  depth  can  be  expected  to  be  of  the  order  of  Ks(Tvap  -  Tj/q"  . 

The  heat  flux  q"g  is  transmitted  forward  in  the  gas  phase  to  heat  the  solid. 

If  ug  is  the  opposing  gas  flow  velocity,  the  forward  transmission  is  described 

by  a  Stanton  number  similarity  mechanism  so  that  q"g  «  ug  Cpg(Tf  *  Tvap). 

Its  delivery  to  the  solid  in  the  vicinity  of  the  leading  edge  Is  described  by 
the  essential  relation  q"g  -  K  (Tf  “  Tvap^V  Substitution  of  these 
relations  into  the  demand/supply  ratio  governing  the  spread  process  leads  to 


»sCpsV  <Tvap  - 

pgcpgKgug  ^Tf  "  Tvap^ 
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This  relation  is  identical  to  de  Ri s *  thick  body  result  in  the  absence  of 
radi ation. 


8.3  Upward  Fire  Spread  on  Vertical  Surfaces;  This  problem  of  practical 

importance  Is  characterized  primarily  by  the  fact  that  the  gravitationally 
induced  hot  gas  flow  is  in  the  same  direction  as  the  flame  spread.  As  a 
consequence  of  the  enhanced  heat  transfer  to  the  upstream  fuel  which  is  yet  to 
be  involved  in  the  upward  spreading  flame,  the  flame  spread  rate  is  higher. 

It  is  evident  from  practical  observation  that  as  the  fire  grows  progressi vely 
larger,  its  growth  rate  also  progressively  increases  as  prompted  by  the 
progressively  greater  heating.  Estimations  of  this  accelerative  growth  places 
it  between  an  exponential  (exp  t)  to  a  power  law  (tn)  type  time-functions  for 
the  flame  spread  rate. 

While  the  mechanism  of  upward  spread  rate  is  qualitatively  known  from 
practical  observations,  the  available  scientific  analysis  of  this  problem  is 
relatively  limited.  Pel  1 o  [19]  reports  an  analysis  of  upward  laminar  fire 
spread  while  Orloff  et  al  [20]  report  on  certain  turbulent  fire  spread 
studies.  Pello's  analysis  points  out  that  while  the  thermodynamic  aspects  of 
the  problem  are  describable  by  the  Stefan  number  CpgTg^L  and  Spalding's 
B-number  [f  Y0a>  hc  +  Cpg(Ts  -TM)]/l,  the  fluid  mechanical  aspects  are  governed 
by  the  parameters 
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where  L  is  (J/kg)  is  the  latent  heat  of  pyrolysis  of  the  solid,  f  is 

stoichiometric  fuel/oxygen  mass  ratio,  Y0w  is  oxygen  mass  fraction  in  the 

ambient  air,  Ts  Is  vaporizing  surface  temperature  (same  as  Tvap  in  de  Ris' 

results),  g  is  acceleration  due  to  gravity,  xjj  is  the  length  of  the  fuel 
surface  between  the  pyrolysis  and  burn-out  stations,  Is  ambient  gas 

kinematic  viscosity,  u  is  rate  of  spread  of  the  pyrolysis  front  and  j>  is 

distance  from  the  origin  fixed  at  the  original  ignition  point. 

Sibulkin  and  Kim  [21]  peformed  an  analysis  in  which  a  parameter 
<t>  a  Vp/ri^hc  is  defined  with  q'F  representing  the  energy  transfer  rate  per 

unit  width  of  the  upstream  uninvolved  fuel  bed  and  m'^  is  the  pyrolyzate  mass 
flow  rate  per  unit  width.  Gravity  enters  into  both  these  quantities.  If 
exceeds  a  critical  value  of  Cps(Tp  -  Tj/hc,  the  flame  spread  rate  is 
unsteady,  accelerative.  In  the  steady  spread  regime,  the  following  relation 
is  deduced. 

PsCRsKA(-1.P.  Lll.-V  =  — 1— r  (8.12 

q'(xp)q"(xp)  (4/3)J 
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j  -  1  for  laminar  spread  and  j  =  0  for  turbulent  spread.  q'c  =  m'fhc  is  the 
combustive  heat  release  rate  (W/m)  per  unit  width  of  the  fuel  plate.  (m'f  is 
the  fuel  production  rate  per  unit  width.)  q"(xp)  is  the  energy  feedback  rate 
(W/m2)  at  the  pyrolysis  front  x  a  xp.  This  flux  is  relatable  to  the  pyrolysis 
rate  m"  at  x  =  xp  through  the  latent  heat  of  vaporization  L  as  q“(xp)  = 
m"(xp)  xL.  The  nondimensional  flame  spread  velocity  parameter  on  the 
left-hand  side  of  the  foregoing  relation  bears  a  remarkably  close  similarity 
with  that  derived  in  the  proceeding  section  for  horizontal  wind-opposed  flame 
spread. 

It  is  thus  evident  that  although  the  focus  of  the  existing  works  on  the 
upward  flame  spread  problems  is  predominantly  on  the  involved  mechanisms, 
certain  scaling  relations  can  be  deduced  from  their  basis.  It  must  be  noted 
however,  that  when  large  scale  fires  are  of  concern,  thermal  radiation  effects 
come  into  play  profoundly  and  can  not  easily  be  incorporated  in  the  existing 
framework  of  knowledge. 


8.4  Fire  Spread  In  a  Forest:  Many  attempts  have  been  made  to  model  the 

problem  of  fire  spread  in  a  forest.  Some  salient  steps  of  progress  will  be 
discussed  here  to  adapt  their  results  for  the  subject  scaling  endeavors. 

Vogel  and  F.  A.  Williams  [22]  experimented  with  fire  spread  across  a 
linear  discrete  array  of  vertically  oriented  match  sticks  In  which  the  scale 
is  so  small  as  to  make  convective  heating  more  important  than  radiative 
heating.  If  the  fuel  elements  are  thick  and  their  spacing  is  small ,  the  flame 
jumps  from  one  element  to  the  adjacent  one  In  a  time  much  shorter  than  the 
characteristic  conduction  time  of  the  solid.  Thus  solving  the  element 
conduction  equation  with  a  convective  boundary  condition,  the  time  to  ignition 
of  adjacent  stick  is  obtained.  Denoting  the  stick  spacing  by  s  and  ignition 
time  by  t-j ,  the  average  fire  spread  rate  u  *  s/t^.  The  result  of  Vogel  and 
Williams  can  then  be  recast  into  the  form 


2 

(8.1: 

where  y*  is  the  flame  standoff  distance  on  the  burning  element  (so  that  (s 
-y*)  is  the  distance  between  the  flame  and  the  next  ignited  stick),  is  the 
stick  thermal  diffusivity  (05  3  Ks/ps^ps)»  K  thermal  conductivity 
(subscript  g  for  gas  and  s  for  the  fuel  solid),  Tg  is  hot  combustion  gas 
temperature,  T0  is  initial  temperature,  and  Ti  is  temperature  of  the  solid 
surface  at  the  instant  of  inception  of  its  flaming.  The  quantity  y*  depends 
upon  the  stick  height  1  as  y*  «  ln  where  n  =  1/4  for  laminar  flames  and  n  *  0 
for  turbulent  flames. 

If  the  stick  burnout  time  t^  is  smaller  than  the  time  to  ignition  ti  of 
the  adjacent  stick,  Vogel  and  Williams  expect  that  the  fire  would  cease  to 
spread.  It  is  observed  experimentally  that  a  particular  stick  ignites  at  its 
top  and  the  flame  spreads  downward  at  a  speed  of  r  cm/sec.  The  stick  is 
generally  noted  to  be  consumed  by  the  time  the  flame  reaches  near  the  bottom. 
Thus,  Vogel  draws  tfo  «  i/r.  Thus  the  fire  spread  is  expected  to  cease  if 
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The  rate  of  downward  spread  may  be  obtained  by  the  results  discussed  in  the 


preceeding  section.  That  is 
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where  ig  is  the  gas  phase  conduction  characteristic  length  which  is  of  the 
same  order  of  magnitude  as  the  combustible  mixture  quenching  distance. 
Substitution  gives  the  condition  for  ceased  fire  spread  as 
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The  temperature  ratio  appearing  on  the  right-hand  side  of  this  criterion 
relation  is  expected  to  be  larger  than  unity  but  only  slightly.  Recalling 
that  y*  «  4n  and  that  a  combustible  mixture  of  high  energetic  Intensity  and 
reactivity  will'  exhibit  a  smaller  quenching  distance  £g,  it  is  clear  from  the 
given  relation  that  short  fuel  elements  spaced  further  will  fail  to  support  a 
spreading  fire. 

Emmons  and  Shen  [23]  studied  the  horizontal  fire  spread  over  a  fuel  bed 
made  of  paper  strips  of  height  H  placed  parallel  to  one  another  at  a  spacing 
of  S.  The  flame  spread  rate  u  along  the  strips  has  been  nondimensional ized 
with  the  paper  thickness  t,  paper  density  p  and  conductivity  K. 
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This  has  been  shown  to  depend  upon  the  ratio  of  strip  height  to  spacing  ratio 
H/S.  When  H/S  <  2,  laminar  flames  are  observed  (with  their  height  If  «  upiH) 
and  when  H/S  >  2,  turbulent  flames  are  noticed  (Lf  *  (uH/g 1 /2S) 2' 3,  where  g  is 
acceleration  due  to  gravity.  Furthermore,  the  flames  on  adjacent  paper  strips 
would  merge  if  the  spacing  is  small,  i.e.,  if  S  <  0.4  +  0.26  H  with  S  and  H  in 
cm. 

While  Vogel's  study  accounts  for  only  convective  heat  transfer  from  the 
hot  flame  gases  to  the  element  yet  to  be  ignited,  Shen's  study  acounts  for 
radiant  heat  transfer  from  the  embers,  convective  heat  transfer  from  the  flame 
gases  and  radiative  heating  from  the  flames  with  an  assumed  flame  height 
correlation. 
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Ataliah  [24]  proposed  a  forest  fire  model  involving  radiation  to  arrive 
at  the  following  result 


where  A,  3,  C  and  D  are  scaled  variables  defined  below. 

A  3  —  fire  spread  velocity  (8*1 
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radiant  heating. 
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H  here  is  the  height  of  the  flar'. 


a  and  Ks  are  respectively  the  thermal 


diffusivity  and  conductivity  of  the  fuel  bed  whose  thickness  is  x.  h  is 
convective  heat  loss  coefficient  between  the  preheating  fuel  and  the 
surroundings.  es  and  £f  are  the  emissivities  of  the  fuel  bed  and  the  flame. 
a  is  Stefan-Bol tzmann  Constant.  T»,  Tp  and  Tf  are  temperatures  respectively 
of  the  surroundings,  the  pyrolyzing  fuel  and  the  flame. 

The  flame  height  H  in  this  scaling  argument  is  yet  an  unsealed  dependent 
parameter.  In  order  to  obtain  this  scaling,  we  resort  to  the  available  flame 
height  analyses  as  following. 

The  inlet  momentum  of  the  fuel  supply  in  buoyant  flames  is  extremely 
small;  the  momentum  produced  by  buoyancy  becomes  the  all  important  issue  in 
determining  the  size  of  a  flame  and  the  wind-flow  induced  by  it.  Thomas 
[25-27]  studied  the  problem  by  applying  dimensional  analysis  methods  to 
correlate  the  experimentally  measured  heights  of  turbulent  stationary  fires. 

Thomas  argues  the  height  of  a  turbulent  buoyant  fire  is  expectable  to 
depend  upon  the  volumetric  fuel  feed  rate  V,  for  it  is  a  measure  of  not  only 


the  total  quantity  of  fuel  to  flow  through  the  plume  but  also  of  the  amount  of 

air  arawn  Into  combustion.  Since  the  velocity  Induced  Ug  in  a  gravitational 

field  of  acceleration  g  by  a  fire  of  base  diameter  d0  has  to  be  of  the  order 

of  /(g  do)  (so  that  the  Froude  number  u2/gd0  is  preserved),  and  since 
•  2  9 

ug  ~  V/do  *  flame  height  H  measured  in  units  of  dQ  is  expected  to  depend 

upon  the  nondimenslonal  ratio  of  either  u^/gdQ  or  The  source  of 

buoyancy  in  a  fire  is  the  lower  density  of  the  fire  gases  compared  to  the 
ambient  air.  This  density  difference  is  mainly  a  consequence  of  the  high  flame 
temperature.  To  account  for  the  larger  buoyancy  produced  by  hotter  fires, 
therefore,  the  dimensionless  temperature  difference  &/Tf  -Tj  has  to  enter 
the  correlation.  0V  is  the  volumetric  expansion  coefficient  (1/K)  of  the  fire 
gases.  Since  body  forces  are  associated  with  the  gravitational  constant,  this 
expansion  term  enters  into  the  denominator  of  the  Froude  number.  Thus 
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The  velocity  utip  attained  by  the  buoyantly  accelerated  flame  gases  at  the 
flame  tip  has  to  be  a  function  of  both  the  flame  height  and  the  modified 
Froude  ratio. 


''tip  a  f  (H  _  V2  \ 
W  2  Vo  ’  g  dQ5  9^  at) 


(8.24) 


The  fuel  input  ^ate  is  sometimes  known  in  terms  of  mass  rate  m0  (kg/s)  which 
is  related  to  the  volume  flow  rate  V  as  V  =  mQ/p0*  Thus  the  Froude  number  i-« 
the  preceeding  two  equations  can  take  the  following  form 
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Since  the  energy  release  rate  Q  in  the  fire  can  be  related  to  the  fuel  feed 
rate  m0  through  the  heat  of  combustion  hc  as  ^  =  ^0hc»  yet  another  way  of 
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expressing  the  Froude  number  exists. 


~r~r- 

po  hcg 


0vAT 


(8.26) 


Throughout  these  dimensional  deductions,  consideration  is  not  given  to  the 
fact  that  the  fire  fluid  elements  at  different  heights  in  the  flame  experience 
different  levels  of  buoyancy.  Steward  [28,29]  analyzed  the  distributed 
buoyancy  problem  to  arrive  at  the  following  prediction  for  flame  height. 
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where  k  is  jet  entrainment  coefficient  («  0.10),  r  is  the  stoichiometric 
air/fuel  mass  ratio,  and  w  s  T00/(T  -  Tj  *  p/(p„"  p)  is  the  inverse  of 
volumetric  expansion  due  to  combustion,  T  and  p  being  temperature  and  density 
of  the  flame  gases.  In  addition  to  being  a  predictive  relation  for  the 
heights  of  stoichlometrically  burning  buoyant  turbulent  fires,  this  equation 
of  Steward  retains  the  basic  form  of  the  Froude  number  given  by  the  earlier 
equation  with  only  the  expansion  term  spelled  out  by  the  function  of  r  and 
■j  square-bracketted  above. 

The  available  flame  height  data  when  compared  to  the  above  prediction 
indicates  that  while  the  predicted  1 /5th  power  dependency  is  satisfactorily 
verified,  the  mixtures  existing  in  real  fires  are  clearly  not  stoichiometric. 
Steward  estimates  that  an  assumed  400%  excess  air  would  make  his  predictions 
quite  close  to  the  experimental  data.  Experience  with  natural  fires  indicates 
that  this  excess  air  figure  is  quite  expectable. 

These  stationary  fire  height  correlations  might  be  adaptable  to 
propagatiro  fires  if  the  fire-seat  diameter  d0  and  the  fuel  inlet  velocity 
(Q/h^ t  dQ  o  )  can  be  suitably  redefined.  Oenote  the  mass  of  fuel  per  unit 
area  of  the  forest  floor,  that  will  be  burned  by  the  fire  sweeping  at  a  speed 
of  u,  by  the  symbol  m".  Let  m"  he  the  average  burning  rate  per  unit  area. 
While  m"(kg/m2)  is  a  prescribable  fuel-load  characteristic  of  the  forest, 
m"(kg/m'sec)  will  depend  upon  the  thermochemistry  and  combustibility  of  the 
forest  fuel.  Then  the  width  w,  of  the  burning  zone  of  the  spreading  fire, 
will  be  given  by  w  *  u  m"/m".  This  w  is  generally  expectable  to  be  nearly 


36 


linearly  proportional  to  the  forest  fuel  bed 
that  d0  In  Steward's  correlation  formula  can 
the  fuel  feed  rate  (Q/p  nd'lh  )  has  obviously 

U  U  U 

doing  this,  the  following  set  of  forest  fire 
realized  from  Atallah's  work. 


height  t.  It  appears  reasonable 
be  replaced  hy  w.  Furthermore, 
to  be  replaced  by  (m"/p  ).  Upon 
spread  scaling  parameters  are 


urn’ 


t  m 


fire  spread  velocity 


•s-  _  Kt 
8  - 


convective  loss 


oe  (T  4  -  T  4) 
■q  =  S  p  °o 

h(T_  -  T  ) 

p  <*>' 


radiative  loss 


D  = 


"SMTf4  - T.4> 

*<vu 


radiative  heating 


E  *  1.49 


(  1  m11  ^  1 

\  1  /  m  \  1 

r(rT./To'  + 

|7(V  * 

(l-„)5  J 

Plume  dynamics 


F  = 


2-.. 

x  m 

a  m,r 
S 


fuel  bed  structure 


(8.28) 

(8.29) 

(8.30) 

(8.31) 

(8.32) 

(8.33) 


Emmons  [30]  presented  an  illuminating  analysis  of  a  radiatively  driven 
forest  fire.  A  fine  spread  rate  parameter  R  and  a  burning  characteristic  P 
are  defined  as  following. 
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The  symbols  u,  m"  and  m",  as  before,  respectively  stand  for  the  flame  spread 
rate,  fuel  loading  per  unit  forest  area  and  average  burning  rate  per  unit 
area.  These  quantities,  as  before,  determine  the  width  of  the  burning  zone, 
w  =  um"/m".  The  factor  a  Is  a  descriptor  of  the  structure  of  the  fuel  bed, 
typically  a  surface  to  volume  ratio,  and  has  units  of  1/meter.  Cp  is  the  fuel 
specific  heat;  and  T0  are  respectively  the  ignition  and  initial 
temperature  of  the  fuel  solid.  Qb  is  the  radiant  emission  Intensity  of  the 
burning  zone.  If  eb  is  the  burning  zone  emissivity  and  Tb  is  its  temperature, 
Qb  is  given  by  (a  ebTb4t)  where  a  Is  the  Stefan-Boltzmann  Constant  and  r  is 
the  forest  fuel-bed  height. 

Emmons  shows  that  only  P  <  1  presents  a  spreading  flame.  As  P  decreases 
from  unity  to  zero,  the  fire  spread  rate  R  and  the  fire  zone  thickness 
w  ■  (R/a)  increase  from  zero  to  a  large  value.  When  R  is  large,  R  «  1/P. 

Addressing  the  issue  of  a  firebreak  of  width  g  in  a  forest  of  height  t, 
it  is  deduced  that  a  radiatively  driven  forest  fire  spread  is  curtailed  if 
g/x  '  R.  Furthermore,  if  w^j  is  the  width  of  the  initial  ignition,  the  fire 
would  grow  to  sustain  and  propagate  only  if  a  wi  >  -  nn(l  -  P).  That  is,  a 
small  P,  indicative  of  high  combustibility  of  the  forest  fuel,  will  support  a 
sustained  fire  even  with  a  ’•'eUti  vely  narrow  ignition  strip.  In  contrast,  if 
P  is  larger  (near  unity),  the  fuel  is  poor  in  combustibility  and  the  ignition 
width  wi  may  have  to  be  larger  than  the  steady-spreading  fire  width  w  In  order 
to  produce  a  sustained  fire. 

While  Emmons  dea.s  with  a  number  of  other  interesting  aspects  of  the  fire 
development  in  a  forest,  it  Is  sufficient  for  us  to  note  the  following  scaling 
parameters  suggested  by  this  work  for  the  radiatively  controlled  forest  fire. 

R  s  a  m"u  (8.31 
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G  =  g/t 
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The  important  effect  of  wind  on  flames  is  not  explicitly  discussed  in  any 
of  the  above-mentioned  Investigations.  A  subtle  exception  to  this  statement 
may  be  Steward's  work  which  shows  that  the  upward  buoyancy  and  the  horizontal 
induced  wind  inertia  conspire  to  cause  a  narrowing  of  the  flame  at  a  small 
height  above  the  fuel  source.  But  more  relevantly,  Welker,  Pipkin  and 
Sliepcevich  [31]  have  studied  the  bending  of  pool  fires  in  a  cross-wind. 

Hamada  [32]  Investigated  the  bending  of  building  fires  in  cross-wind.  Both  of 
these  works  indicate  that  the  Froude  and  Reynolds  numbers  are  the  primary 
criteria  governing  the  flame  tilt  angle.  Welker's  momentum  balance  on  the 
flame  leads  to  the  result 


(8.40) 


where  e  is  tilt  from  vertical;  K  Is  a  constant  (near  about  60);  f  is  a  flame 

shape  factor  (upright  cylinder  f  3  tt/4),  slab  f  3  1/2,  cone  f  =  it/6,  pyramid 

f  =  2/3,  and  paral lelepided  f  3  1);  pf  and  Pa  are  the  densities  respectively 

of  the  flame  gases  and  ambient  air;  v  is  wind  speed;  g  is  gravity;  d  is  fire 

base  diameter  and  u  is  ambient  air  dynamic  viscosity.  The  group  v^/gd  is 
Froude  number  while  p  vd/u,  is  Reynolds  number.  The  indices  n  and  m  are  equal 

a  a 

to  1/2.  Hamada's  experiments  lead  to  a  similar  result  with  the  mere 
difference  that  n  »  1  and  m  -  0. 

Welker  and  Sliepcevich  [33]  reported  that  the  temperature  in  the 
convective  plume  tilted  by  the  wind  can  be  uniquely  correlated  in  the 
framework  of  the  following  parameters  for  wide  line  fires. 


(8.41) 


(8.42) 
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where  T  Is  temperature  In  the  plume  at  a  distance  x  downstream  from  the  center 
of  the  pool,  Ta  is  ambient  temperature  and  O'  is  total  heat  output  rate  of  the 
fire  per  unit  length,  w  is  width  of  the  fuel  pool  channel,  z  is  height  from 
the  floor.  Attention  is  drawn  to  the  manner  in  which  the  heat  release  rate 
defines  the  characteristic  velocity  in  this  context.  A  recent  study  by  Mudan 
and  Croce  [34]  also  confirms  the  importance  of  the  Froude  number  in 
determining  the  characteristics  of  wind-blown  pool  fires. 

Finally,  our  discussion  of  forest  fire  modeling  will  be  Incomplete  if 
mention  is  not  made  of  Rothermel's  work  [35],  With  fuel,  weather  and 
topographic  parameters  measurable  in  the  field  as  inputs,  the  model  is  capable 
of  predicting  both  the  forest  fire  spread  rate  and  the  fire  intensity. 

Keeping  consistent  with  our  earlier  notation  but  retaining  the  British  units 
of  the  original  work,  Rothermel's  equation  for  fire  spread  rate  u  (ft/min)  is: 

_ (8.43) 

m"  5(1  +  ^  +  <(>s) 

where  pb  is  fuel  bulk  density  (lbm/ft3);  e  is  an  effective  heating  number 
(e  =  exp(-138/<j)  where  a  is  fuel  particle  surface  area  to  volume  ratio, 

1/ft);  Qi  is  energy  of  preignition  (Btu/lbm)  (Qi  -  Cp ( T-j  -  T0))  In  the 
perspective  of  our  earlier  discussions  where  Cp  is  the  fuel  specific  heat, 
is  its  Ignition  temperature,  and  T0  is  its  initial  temperature  and  0-j  =  250  + 
1,116  Mf  for  forest  fuels  of  moisture  content  Mf  lbm  of  moisture/lbm  of 
oven-dry  wood);  m"  reaction  intensity  (Btu/ft^min)  (a  function  of  optimum 
chemical  reaction  velocity,  net  fuel  loading,  fuel  particle  heat  content, 
damping  of  the  reaction  due  to  moisture,  and  damping  of  the  reaction  due  to 
mineral  content);  £  is  propagation  flux  ratio  (a  function  of  the  fuel  particle 
surface  to  volume  ratio  and  of  the  fuel  packing  ratio,  3);  <jw  is  wind  factor 
(a  function  of  the  wind  speed,  a,  and  3);  and  ^  is  the  slope  factor  (a 
function  of  6  and  the  slope  of  the  terrain).  Note  that  e,  £,  ^  and  ^  are 
all  nondimensi oral  factors. 

Rothermel's  work  has  enjoyed  remarkable  success  in  forest  fire  control. 

It  is  interesting  to  note  that  the  product  of  Emmons  parameters  P  and  R, 

PR  s  in"  Cp ( T-j  -  T0)u/Qb,  shares  a  remarkably  common  conceptual  basis  with 
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Rothermel's  fire  spread  rate  parameter  in  the  form  presented  above.  The  main 
difference  is  that  Rothermel's  model  goes  into  much  greater  detail  of  the 
burning  zone  heat  flux  than  does  Emmons'  model. 

8.5  Fire  Growth  In  Compartments t  While  the  previous  section  has  dwelt 
much  on  fire  spread  outside  buildings,  we  deal  in  this  section  with  the 
plausibility  of  scaling  the  fire  growth  process  within  a  room.  A  recent  paper 
by  Kanury  [4]  forms  the  basis  of  this  discussion  although  literature  is 
plentiful  of  Investigations  on  this  issue. 

As  the  Initial  fire  burns  under  conditions  of  essentally  unrestrained  air 
supply,  its  hot  gases  and  flames  convect  and  radiate  energy  to  the  enclosure 
walls  In  a  manner  affected  by  the  fire  location/size  and  the  room/ventil ation 
details.  Conductive  heating  of  the  linings  will  lead  to  pyrolysis, 
vaporization  or  charring  to  produce  combustible  vapors  which  then  convecti vely 
mix  with  air  and  ignite  to  enhance  the  flame  volume.  The  conductively  heated 
walls  may  emit  radiation  to  make  the  room  a  sort  of  an  energy  augmentor,  an 
Introspective  hohlhrum.  Accumulation  of  smoke  near  the  ceiling,  however,  may 
soon  attenuate  the  surface  radiant  heat  exchange  and  introduce  soot  layer 
radiation  into  the  subsequent  heating  process.  All  during  this  transient 
process;  either,  the  initial  fire  is  continuously  and  acceleratively  growing 
towards  the  state  of  a  physico-chemical  catastrophe  know  as  flashover;  or,  if 
the  initial  fire  is  externally  controlled  to  maintain  a  constant 
size/intensity  such  as  could  be  done  with  a  gas  burner,  the  room  wall  linings 
are  gradually  heated  to  arrive  at  a  critical  thermochemical  state  at  which  the 
various  surfaces  tend  to  be  simultaneously  ignited. 

The  variables  of  influence  identifiable  by  this  description  exhibit 
nonlinear  interactions  among  the  various  time-dependent  fire  processes.  The 
room  geometry,  for  instance,  mutilates  the  thermal  radiative  exchange  of 
energy  in  the  room  as  well  as  alters  the  convective  flows  of  air  and  hot  gases 
and  the  associated  mixing  and  heating  phenomena.  As  a  result,  the  ignition 
(and/or  fire  spread)  will  be  influenced  by  the  room  geometry. 

The  ventilation  intensity  will  obviously  affect  the  air  supply  rate  to 
the  fire,  convection  and  mixing,  accumulation  of  smoke  in  the  ceiling  layer 
and  hence  the  overall  thermal  radiation  process  itself.  Herein  arises  the 
important  connection  between  the  window  geometry  (which  determines  the  flow  of 
gases)  and  the  lintel  height  (which  determines  the  thickness  of  the  hot 


ceiling  gas  layer).  The  energy  balance  responsible  for  continued  fire-growth 
process  (or  the  ignition  process  of  various  parts  of  the  room)  is  thus 
tampered  by  the  ventilation  character! sties. 

The  nature  of  the  Initial  fire  and  its  proximity  to  the  lined  walls  shall 
result  In  convective  and  radiative  heating  of  the  adjacent  surfaces,  thus 
preparing  them  towards  the  thermochemical  state  of  ignition.  In  this  entire 
picture  of  heat  exchange,  the  wall  linings  play  not  only  the  active  role  of 
confining  the  gas  flows  and  absorbing/reflecting/emitting  thermal  radiation 
but  also  a  passive  role  of  being  heated  by  transient  conduction  progressively 
towards  a  thermochemical  state  of  receiving  the  growing  fire  or  of  spontaneous 
ignition. 

This  qualitative  discussion  leads  one  to  believe  that  the  flashover 
problem  entails  two  specific  questions:  (a)  Is  flashover  possible  under  a 
given  set  of  conditions  of  the  room  geometry,  ventilation,  initial  fire  and 
wall-lining?  (b)  If  it  is  possible,  how  does  the  time  to  flashover  depend 
upon  these  prescribed  conditions?  The  answers  to  these  questions  require  that 
the  experimental  fire  growth  data  be  represented  on  the  basis  of  a  dimensional 
analysis  of  the  problem  of  room  fire  growth  and  the  associated  transient 
heating  of  the  room.  It  is  apparent  that  the  variables  expected  to  influence 
the  growth  of  a  room  fire  to  flashover  pertain  to  the  following  five  groups: 

a.  room  geometry  —  volume,  depth  and  height,  lintel  height,  etc; 

b.  ventilation  —  door/window  location  and  dimensions; 

c.  fuels  --  furnishings  and  wall  linings; 

d.  initial  fire  —  size,  location,  burning  rate,  etc;  and 

e.  heating  of  walls  and  ceiling  by  conduction,  convection  and 
radi ation. 

If  the  variables  involved  are  thus  identified,  five  independent 
dimensionless  parameters  are  expected  to  evolve  pertaining  to  the  five 
families  of  variables  listed  above.  (For  the  convenience  of  this  discussion, 
these  parameters  will  be  denoted  by  A,  B,  C,  D,  and  E.)  Typically  then  one 
can  represent  the  time-wise  development  of  the  enclosure  fire  by 
systematically  altering  one  of  the  parameters  (say,  the  initial  fire 
parameter)  while  the  other  four  parameters  are  kept  fixed.  A  critical  value 
of  the  altered  parameter  is  expected  to  exist  representing  the  occurrence  of 
flashover.  Beyond  the  critical  value  of  the  parameter,  the  time  to  flashover 


42 


will  depend  on  the  magnitude  of  the  parameter. 

Kanury  [4]  deduces  these  scaling  parameters  from  the  laws  of  conservation 
of  mass,  energy  and  momentum.  Whereas  mass  and  energy  equations  dominate  in 
the  deduction,  the  momentum  equation  comes  into  play  only  indirectly  in 
determining  the  various  mass  fluxes  involved.  The  species  equations  are  not 
considered;  in  the  estimation  of  thermal  radiation,  smoke  is  taken  as  a 
passive  variable  its  generation  rate  being  considered  as  a  property  of  the 
fuel  ensemble. 

The  mass  conservation  requires  under  quasi-steady  state  that  the  mass  of 
air  (ma)  and  fuel  (m  )  entering  into  the  control  volume  around  the  room  be 
equal  to  the  mass  of  combustion  products  leaving  the  room  (mg).  While  mv  will 

be  determined  by  the  initial  fire  burning  characteristics,  and  m„  mutually 

a  g 

adjust  to  share  the  window  opening  for  flow  according  to  the  well-known 
Fujita-Kawagoe  equations  In  which  the  window  height  H0  and  width  W0,  densities 
of  the  combustion  product  gas  and  ambient  gas,  gravitational  constant,  and  the 
window  orifice  coefficient  appear  as  prescribed  parameters.  Coupling  these 
relations  with  the  relation  describing  the  entrainment  of  induced  air  into  the 
fire  plume  to  mix  and  react  with  the  fuel,  one  obtains  a  fully  posed  problem 
of  flow  in  the  room.  The  solution  of  this  problem  yields  the  scaled  airflow 
rate  *  as  a  function  of  the  (prescribed)  scaled  fuel  input  rate  <f>  and  of  the 
window  aspect  ratio  i|>:  4=  <»(  ^ ,  ip )  with  the  following  definitions. 


*  3/2 

» "  'VkVo 

(8.44) 

*  ■  'VkVo3/? 

(8.45) 
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k  is  a  constant  («  1.6  kg/m5/2s)  composed  of  the  product  gas  and  ambient  air 
densities,  the  gravitational  constant  and  the  window  orifice  coefficient.  The 
location  of  the  initial  fire  within  the  room  is  an  additional  issue  which 
exerts  a  minor  influence  on  the  dependency  of  $  on  $  and 

The  energy  released  due  to  (total  or  partial)  combustion  of  the  fuel  mass 
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supply  is  used  in  part  to  raise  the  enthalpy  of  gas  flow  and  to  heat  up  the 

enclosure  walls  while  the  remainder  is  radiated  out  through  the  window.  (In 

the  preflashover  stage  of  fire  growth,  all  the  combustion  may  be  presumed  to 
be  contained  within  the  enclosure  so  that  flames  do  not  shoot  out  of  the 
opening.)  The  fuel  mass  flow  rate  mv  is  generally  so  small  that  mg  «*ma.  One 
can  ignore  the  radiative  loss  through  the  window  for  the  introspective 
consideration  of  fire  growth  to  flashover.  The  energy  received  by  the 
enclosure  surfaces  can  be  accounted  to  heat  the  (thick)  linings  by  transient 

conduction.  Gas  phase  phenomena  could  be  construed  to  follow  a  quasi-steady 

behavior. 

A  synthesis  of  these  ideas  leads  one  to  arrive  at  a  dependency  of  the 
enclosure  surface  temperature  9w  on  time  t»  air  flow  rate  %  fuel  feed  rate 
ceiling  soffits  skirt  area  y,  and  the  wall-heating  characteristic  n.  Since 
$  =  *((j>,ip),  thus  6w  -  0w(t;  4,  ij>,  y,  n).  The  temperature  of  the  hot  gases 
near  the  ceiling  0g  also  depends  on  the  same  parameters  in  a  nearly  same 
functionality.  The  scaled  quantities  are  defined  as  follows: 


,(T<  -  TJ 


j  a  g.w 


T  =  ^t/K  p  0 


wpw  pw 


y  s  7T/W0K0 


kC  H 

pg  o 

Subscripts  g  and  w  respectively  represent  the  hot  gas  and  the  heated  surface, 
a  is  for  ambient  air,  o  is  for  the  opening.  While  t  is  time,  T  is 
temperature.  K  is  thermal  conductivity,  k  is  window  flow  constant  mentioned 
earlier.  p  and  Cp  are  density  and  specific  heat  respectively.  W  is  width  and 

H  is  height.  hc  is  heat  of  combustion  of  the  fuel.  A  is  the  area  of  the 
room  surface  in  its  upper  reaches  which  will  become  contacted  by  the  hot  smoke 


layer;  thus,  win  be  the  celling  area  plus  the  product  of  the 
room-wal 1 -perimeter  near  the  ceiling  and  the  smoke  layer  depth,  h*  is  the 
combined  convective  and  radiative  heat  transfer  coefficient  between  the  hot 
gases  and  the  enclosure  surfaces  exchanging  heat  with  these  gases.  The 

radiative  component  of  TT  is  yet  to  be  fully  investigated. 

If  the  initial  (prompting)  fire  characterized  by  $  is  situated  adjacent 
to  a  wall,  the  wall  linings  are  expecxtad  to  get  ignited  and  to  burn  so  as  to 
result  in  a  time-wise  growing  prompting  fire.  The  effect  of  this  will  be  to 
expediate  the  flashover.  An  additional  parameter  is  required  to  account  for 
this  effect.  This  parameter  has  to  be  a  measure  of  the  combustibility  of  the 
lining  material  more  specifically,  of  its  thermochemical  response.  If  h  is 
the  heat  transfer  coefficient  between  the  flame  gases  at  temperature  TM  and 
the  lining  at  temperature  Tw,  q  is  heat  release  rate  of  the  lining  sample 
measured  upon  exposing  to  an  irradiance  of  1  and  f  Is  the  empirical  constant 
ratio  of  area  covered  by  a  flame  to  its  power  (f  ~  0.0?  m2/kW),  then  a  wall 
lining  thermal  response  parameter  is  definable  as 

8  s  f  h(TM  -  Tw)  q/I  (8.51) 

This  would  then  be  an  additional  parameter  of  the  room  fire  flashover  scaling 
problem. 

An  application  of  these  scaling  arguments  to  the  existing  room-fire-test- 
data  indicated  that  the  approach  is  in  correct  direction  although  a  more 
explicit  account  of  thermal  radiation  appears  to  be  necessary  as  also  a  more 
rigorous  consideration  of  the  wall-lining  combustibility. 

8.6  Radiant  Heating  and  Jumping  Fire  Spread:  While  a  clean  separation 

of  the  topics  of  relevance  is  not  completely  possible,  two  important  elements 
of  existing  knowledge  on  fire  jump  across  a  distance  of  separation  between  two 
buildings  will  be  discussed  here. 

First,  the  study  of  Mudan  and  Croce  [34]  shows  that  the  thermal 
radiation  from  trench  pool  fires  decays  with  distance  according  to  the 
inverse-square-law  from  a  flux  exceeding  about  15  kW/m^  near  the  edge  of  the 
fire  to  about  0.5  kW/m2  at  a  distance  of  about  15  pool  (hydraulic)  diameters. 
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Second,  Margaret  Law  [36]  reported  on  safe  separation  distances  between 
buildings  to  prevent  ignition  of  successive  buildings.  Riven  the  Intensity  of 
radiation  emitted  by  a  fire,  and  the  dimensions  and  distribution  of  windows 
and  other  openings  of  the  burning  building.  Law  shows  that  it  is  possible  to 
calculate  the  distance  beyond  which  the  impingent  radiant  intensity  at  a  point 
on  a  vertical  facade  facing  the  burning  building  would  be  lower  than  that 
required  to  causo  ignition  of  the  facade.  The  problem  becomes  one  of  mere 
geometry  such  that  this  critical  distance  is  given  by  the  equality  in 

N 


where  Fn_f  is  the  radiation  view  factor  between  the  nth  window  opening  of  the 
burning  building  which  has  N  total  openings  and  the  target  facade  f.  Ic  is 
the  critical  irradiance  below  which  the  facade  material  can  not  be  ignited 
(Ic  «  12.5  kW/m2  for  most  cellulosic  and  synthetic  organic  materials)  and  I0 
is  the  emission  intensity  of  the  already  burning  fire  (I0  »  75-150  kW/m2  for 
most  building  fires  depending  upon  the  degree  of  ventilation  to  the  burning 
room).  Safe  distance,  then,  is  such  that  z  Fn_f  >  0.075  to  0.15 
approximately.  Radiation  heat  transfer  literature  (for  example,  Siegel  and 
Howell  [37])  abounds  with  dictionaries  of  view  factors  between  surfaces  of 
different  shapes  disposed  in  different  geometries.  Of  specific  interest  in 
the  problem  of  fire  jump  between  buildings  is  the  view  factor  between  a  window 
(usually,  rectangular  or  square)  at  a  given  height  and  a  facing  large  parallel 
surface  of  the  target  facade. 
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SECTION  9 


MASS  FIRE 

As  radiant  ignition  of  a  spectrum  of  combustibles  occurs,  the  fire  would 
tend  to  grow  slowly.  Soon  would  arrive  the  blast  wave  to  extinguish  some  of 
these  primary  fires  and  to  augment  others  and,  most  importantly,  to 
mechanically  disrupt  structures  and  installations.  The  degree  of  this 
mechanical  disruption  depends  upon  the  distance  from  ground-zero,  its 
consequence  being  two-fold:  first,  piles  of  debris  will  be  formed  within 
which  the  primary  fires  surviving  the  blast  would  continue  to  grow;  and 
second,  secondary  fires  will  be  triggered. 

The  combination  of  circumstances  would  aid  the  fire  to  spread  in  a  manner 
elicited  in  the  preceeding  chapter.  As  the  total  fire  becomes  large,  it  is 
termed  a  ‘mass  fire'.  A  mass  fire  is  said  (Alger  [38],  Parker  [39])  to  be  a 
fire  of  such  a  large  magnitude  within  which  individual  fires  interact  and/or 
merge  with  one  another  to  produce  a  buoyant  convective  flow  pattern  intense 
enough  to  be  capable  of  perturbing  the  ambient  atmosphere.  A  mass  fire  which 
moves  with  the  prevailing  wind  is  usually  known  as  a  conflagration  while  one 
which  burns  stationary  is  termed  a  firestorm.  Fire  devils,  high  winds,  and 
other  severe  effects  are  generally  associated  with  a  mass  fire  to  result  in  an 
enhanced  destructiveness  of  the  fire. 

In  this  chapter,  we  examine  the  plausibility  of  scaling  mass  fires, 
especially  the  firestorms.  The  discussion  includes,  in  the  order,  the 
following  topics:  (a)  Dynamics  of  plumes  over  'small'  turbulent  fires; 

(b)  radiative  effects  in  plumes;  (c)  Induced  wind  flow;  (d)  firebrand 
transport;  (e)  fire  and  wind/fire  storms;  (f)  the  combustion  zone  specifics. 
9.1  Dynamics  of  Plumes  over  'Small'  Turbulent  Fires:  Priestley  and  Ball 

[40]  studied  in  1954,  the  topic  of  continuous  convection  from  an  isolated 
source  of  heat.  The  velocity  and  temperature  fields  over  finite  size  fires 
can  be  deduced  from  the  given  results.  Let  a  fire  of  radius  R0  be  the 
plume-source  of  vertical  velocity  w0  and  temperature  T0.  Let  its  virtual 
point  source  lie  at  z  =  -Zq  below  the  actual  fuel  surface.  Let  R  be  the 
radius  of  the  plume  at  any  height  z  where  the  axial  velocity  and  temperature 
(excess  over  the  ambient  temperature  Tj  be  wmax  and  (T  -  Tjmax.  The  radial 
(r)  distributions  are  found  to  be  given  by  w/wmax  *  (T  -  TM)/(T  -  Tjmax  - 
exp(-p2/2R2) .  Denote  the  scaled  velocity  and  temperature  by 


(3gA/2T(Bc  zQ) 


(T  -  TJ 

9  H  2  4  5TIT3  (9*' 

(2T//3gc4zo5) 

where  g  is  gravitational  acceleration,  A  =  R?  wmax(T  -  Tjmax,  a  constant, 
indicative  of  the  source  energy  strength  and  c  =  dR/dz,  found  from  the  theory 
to  be  a  constant  near  0.1  indicative  of  the  nature  of  turbulence  which 
describes  the  entrainment  and  mixing  in  the  plume.  Then  the  axial  velocity 
wmax  and  temperature  excess  9niax  are  obtained  to  vary  with  height  2  according 
to  the  following  relations: 


z„,3 
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When  W0  is  less  than  (2/3) 1/3,  the  hot  gases  will  accelerate  to  a  velocity 
W  a  (21/3/3W2)(l  -  Wq3)-1/6  at  a  height  of  z/z0  -  (3(1  -  W03))  1/2  and  then 
decelerate  as  a  result  of  mixing  with  entrained  fluid.  If  W0  >  (2/3) 1 y 3,  the 
accelerating  stage  will  be  absent  and  the  motion  will  decelerate  through  out. 
At  sufficiently  large  heights,  the  first  terms  of  the  above  equations  dominate 
and  the  solution  approaches 

v»  ■  <vz)1,!  <9-5 
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In  a  fire,  W0  is  insignificantly  small  so  that  '4ax  “  C(z0/z)  -  (z0/z)3l 

which  yields  the  highest  velocity  of  W  *  0.73  at  a  height  of  z/z0  =  /3 
correspondi ng  to  a  z/R0  of  7.3. 


If  the  atmosphere  is  stratified,  the  stratification  parameter  y  is 
defined  as 


Y  = 


(2T 


A2/3gcV 


1/3  j 


1/4 
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When  dT^/dz  is  positive,  stable  stratification  Is  indicated;  otherwise,  it  Is 
unstable  stratification.  A  small  range  in  y  covers  a  very  wide  range  In  the 
source  strength. 

Under  stable  stratification  conditions,  the  hot  gases  first  accelerate 
(nearly  independent  of  y)  and  then  rapidly  decelerate  to  zero  velocity  at  a 
certain  maximum  height  which  is  smaller  for  large  values  of  y.  That  is,  there 
then  exists  a  ceiling  beyond  which  the  column  does  not  penetrate.  The  ceiling 
height  is  approximately  given  by  (zce-ning/z0)  -  1.9/y3/2.  Under  given 
environmental  conditions,  the  ceiling  height  will  depend  on  the  source 
strength  but,  to  a  fair  approximation  independent  of  the  source  size. 

The  abrupt  decrease  in  plume  velocity  as  the  height  approaches  the 
ceiling  height  is  mainly  due  to  a  reversal  in  the  plume  buoyancy  which  occurs 
at  an  altitude  of  (z/z0)  «  (1  +  1.8/y3) 1 /2  or,  for  small  values  of  y,  at 
z  -  0.71  Zceiling* 

Under  unstable  stratification  conditions,  the  gases  will  first 
accelerate,  then  decelerate  and  cool  down  but  finally  pick  up  acceleration  to 
ultimately  reach 
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dz 
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Morton,  Taylor  and  Turner  [41]  in  a  classic  paper  published  in  1956 
solved  the  same  problem  as  Priestley  and  Ball  but  with  a  point  source.  Hefine 
the  constant  of  buoyancy  injected  into  the  plume. 

1  =  -  »W».  <9'9 

where  b  is  a  measure  of  the  plume  thickness  (radial  location  at  which  the 
velocity  decays  to  1/2  of  the  umax),  while  b  is  proportional  to  Priestley's 
R,  (p®  -  p)/p»  for  an  ideal  gas  mixture  is  equal  to  (T  -  TJ/T„.  Finally,  0 
is  related  to  Priestley's  A  by  Q  =  (2gA/TM)(b/R)2.  Also  define  a  mixing 
coefficient  a,  as  the  ratio  of  entrainment  velocity  to  the  local  maximum 
velocity  in  the  plume,  to  be  related  to  Priestley's  c  via  a=  (5c/6 ) ( b/R ) . 


Morton's  solution  is 


These  relations  are  essentially  same  as  the  far-field  solutions  of  Priestley 
and  Bal 1 . 

The  measurements  of  Rouse,  Yi h  and  Humphreys  [42]  indicate  that  in  any 
cross  section,  the  plume  velocity  profile  differs  from  the  temperature  profile 
according  to 

~  ■  exp(-rW) 
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where  8^  is  an  empirical  constant  equal  to  96/71.  Using  these  profiles  in  the 
Integral  equations  of  mass,  momentum  and  energy  conservation,  the  following 
results  are  obtained  for  the  weakly  buoyancy  plume  [43]. 


which  are  again  consistent  with  the  results  of  Morton  and  the  far-field 
profiles  of  Priestley  except  for  slight  differences  in  the  coefficients. 
G.  I.  Taylor  [44]  succlntly  discusses  the  implications  of  these  scale 
rel  ati "ns . 


The  scaling  scheme  for  the  plume  dynamics  thus  seems  to  be  In  hand  with 
b  «  (  jZ)  ,  wfflax  *  (Q/a2z)1/3  and  gApmax/p-  *  (Q2/a4z5) 1 '3.  If  the  fire  were  a 
line  rather  than  a  point  source,  then  these  scale  relations  are:  b  «  (az), 
wmax  *  ( Q 1  / a) 1  / 3  and  gApmax/p„  *  (0‘2/a  2^3)1/3  *here  Q1  Is  buoyancy  source 
strength  per  unit  length  of  the  fire  [45]. 

9.2  Radiative  Effects  In  Plumes:  The  relations  of  the  preceeding 

section  make  scaling  possible  of  the  far-field  velocity  and  temperature 
distributions  In  plumes  over  fires  when:  (1)  thermal  radiation  Is  absent;  and 
(11)  the  flame  and  ambient  temperatures  are  same  In  the  prototype  as  in  the 
scale  model.  Written  explicitly  in  scaling  format,  these  relations  are: 


(9.18) 


(9.19) 


wherein,  b/z  =  6<*  /W2  «  0.085  «  a  constant  is  already  incorporated  into  these 

profiles.  Note  then  [6a  /5/2~]2  -  0.0072  and  [6<*b/5 /T]2  -  0.0097.  T5p 
preexponential  (square-bracketted)  constants  are  about  S.6  and  21.2 
respectively  for  the  velocity  and  temperature.  Remember  also  that  the 
buoyancy  source  Q  is  defined  as  b2umax  g(  Pco  -  p)max/ p«>  whose  units  are  m4/s3 
and  is  related  to  the  power  q  of  the  fire  by  q  =  [us'p  C_T  /(g(  s'  +  l))l0. 

oo  p  00 

Further  note  that  the  entrainment  velocity  ue  at  the  edge  of  the  plume  at  any 
height  z  is  given,  by  definition,  as  ue  =  a  w^ax*  Thus,  ue  scales  according 
to 
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Hottel  [46]  notes  that  modeling  thermal  radiation  in  this  natural 
convection  of  the  fire  plume  can  be  done  at  two  levels  of  detail:  the  first 
is  complete  modeling  accounting  for  the  alteration  of  the  flow  field  by 
thermal  radiation;  the  second  is  partial  modeling  only  of  the  radiation  from 
the  plume  to  the  surroundings  with  unaltered  flow  field. 

In  seeking  complete  modeling,  one  attempts  to  preserve  both  the 
convective  modeling  discussed  above  and  the  following  ratio  of  flame/plume 
radiation  to  convective  heat  transfer. 
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(9.21) 
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a  is  Stefan-Boltzmann  constant  and  *  is  the  plume  gas  absorption 
coefficient  (1/m).  The  quantity  *b  is  proportional  to  the  product  of 
atmospheric  pressure  and  the  height.  If  this  product  were  attempted  to  be 
preserved,  Q  has  to  be  proportional  to  z3/2  in  conflict  with  the  buoyancy 
distribution  modeling  requirement  of  Q  *  zs/2.  Such  a  preservation  also 
tampers  with  the  convective  heat  transfer  rate  term  of  the  denominator  of  the 
above  equation  by  calling  for  Q  «  z3,  which  too  is  in  vivid  contrast  with 
nonradiative  plume  scaling  prer.cri ptionr, .  It  is  thus  concluded  that  complete 
modeling  of  radiation  in  flames/plumes  is  impossible. 

If  one  then  contents  oneself  with  the  partial  modeling  in  the  second 
level  of  detail  mentioned  above,  then  one  can  satisfy  the  requirement  of  0  « 
z5/2  if  the  flame/plume  is  radiatively  opaque. 

Murgai  [47]  investigates  the  influence  of  radiative  transfer  on  the 
structure  of  the  turbulent  convective  plume.  Standard  mass,  momentum  and 
energy  equations  are  the  starting  point  O'  the  analysis  but  a  radiative 
heating  rate  Is  added  to  the  energy  equation.  The  effect  of  radiation  is  to 
destroy  tne  constancy  of  the  buoyancy  content  Q  with  respect  to  height,  a 
nondimensional  radiation  number 
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is  defined.  Pw,  y»  I  and  a*  are  nondimensional  and  are  respectively:  the 
ambient  pressure  suitably  nondimensionalized  with  a  reference  pressure, 
exponent  of  the  adiabatic  law,  a  shape  factor  and  a  constant.  All  the 
dimensional  properties  are  in  the  first  bracketed  ratio:  <  is  the  absorption 
coefficient  (1/m)  of  the  fluid,  a  is  Stefan-Boltzmann  constant  (W/rn^K4),  Tw  is 
ambient  temperature  (K),  b0  is  source  radius  (m),  p*  is  ambient  density 
(kg/m3),  Cp  is  plume  fluid  specific  heat  (J/kg  K),  g  is  gravitational  force 
per  unit  mass  (m/s3),  and  (p*  -p)0  is  source  density  defect  (kg/m3).  The 
first  bracketed  nondimensional  term  may  be  interpreted  as  a  ratio  of  the 


characteristic  radiant  flux  glj^  to  a  characteristic  input  convective  heat 
flux  p0OCpT00[g(  p.  -  p)0/( <2boPa#)3l  >z  times  the  fractional  initial  density 
defect  ( p„  -  p) 0/ Poo*  The  quantity  in  the  square-bracket  of  the  characteristic 
convective  flux  is  a  velocity. 

Murgai  sees  that  when  $  is  very  small,  an  optically  thick  (i.e.,  opaque) 
limiting  case  arises  to  reduce  the  problem  to  a  plume  with  no  Internal 
radiative  exchange.  This  is  entirely  consistent  with  Hottel's  conclusion.  If 
$  is  large,  however,  the  plume  Is  optically  thin  (I.e.,  transparent)  to  alter 
its  development.  Murgai  numerically  integrates  the  conservation  equations  to 
obtain  the  spreading  of  the  plume,  and  its  velocity  and  temperature 
distributions.  It  is  found  that  radiation  does  not  noticeably  change  the 
linearity  of  growth  of  the  plume  radius  with  height.  The  more  Intense  the 
radiation,  (i.e.,  the  larger  the  $),  the  fatter  the  plume.  This  not  due  to 
thermal  dilution  ascribable  to  the  mass  addition  by  entrainment  as  in  a  purely 
convective  plume  but  due  to  radiative  cooling.  The  consequent  loss  of 
buoyancy  is  quite  vivid  even  with  minor  ( 4>  =  0.02)  radiation.  So  also  is  the 
rapidity  with  which  the  upward  velocity  decays  with  height. 

Another  important  paper  on  modeling  fire  plumes,  is  that  of  Morton  f 481 
in  which  the  characteristics  of  strongly  buoyant  plumes  as  well  as  the  effect 
of  radiation  on  plume  dynamics  are  discussed. 

Both  the  intermediate  and  lower  regions  of  the  plume  are  generally  much 
hotter  and  Involve  strong  buoyancy  and  thermal  radiation.  Fntrainment  in  the 
strongly  buoyant  regime  is  expected  to  depend  upon  the  ratio  of  local  mean 
plume  density  p  to  the  ambient  density  p„  as  well  as  the  lateral  momentum, 
pw^.  Thus,  remembering  that  the  entrainment  velocity  u  's  equal  to  a  •  w,  a  = 
<*( p/ 0»» ) .  Furthermore,  as  the  plume  cools  down  to  the  weakly  buoyant 
state,  a  +  cvrf,  the  standard  entrainment  constant  for  weak  plumes. 

Additionally,  the  measurements  of  Ricou  and  Spalding  suggest  a*  pW2.  These 
fragments  of  information  led  Morton  to  dimensionally  reason  that  the  most 
viable  as  well  as  convenient  relation  is  a  =  (p/pj1/2 cw*  This  is  consistent 
with  the  common  observation  that  strongly  buoyant  plumes  entrain  significantly 
less  than  the  weak  ones;  and,  as  a  result,  strong  plumes  spread  less  apidly 
than  the  weak. 

Near-field  plume  dynamics  involve  several  closely  coupled  phenomena.  For 
example,  the  reduced  entrainment  may  significantly  influence  the  combustion 
process  itself  which  in  turn  effects  the  rate  of  energy  and  hence  buoyancy 


Injection  Into  the  plume.  (The  fuel  pyrolysis  and  combustion  near  the  source 
are  extremely  complex  to  be  Incorporated  into  a  realistic  scaling  scheme;  see 
the  work  of  Nielsen  and  Tao  [49].)  Radiation  is  Important  <n  a  number  of 
ways.  When  the  absorption  coefficient  <  is  small,  (i.e.,  the  radiation  mean 
free  path  l/«  Is  large  compared  to  the  source  si2e),  the  plume  is  transparent. 
Otherwise,  it  Is  opaque;  then  the  vertical  flux  of  radiation  F  within  the 
opaque  plume  core  Is  related  to  local  'uniform'  temperature  T  through 
4na  dTVdx  =  -  3*  F.  Morton  nondimensionalizes  the  problem  to  realize  the 
following  three  scaling  ratios. 

A  =  8tt  /3k  bQ 

c  *  ra«T04V2%  ’-Vo 

c  5  2  Vo  /sb0 

where  aW  is  the  weak  plume  entrainment  coefficient  (»  0.1),  <  is  absorption 
coefficient,  8  is  coefficient  of  gas  expansion,  and  subscript  o  pertains  to 
the  source  conditions.  A  is  a  measure  of  the  radiation  mean  free  path  length 
to  the  fire  base  radius.  For  opaque  flames,  A  lies  in  the  range  of  10-2  to 
10”4,  so  small  as  to  permit  ignoring  the  vertical  diffusion  of  radiation.  0 
is  a  measure  of  the  radiative  loss  to  the  surroundings  relative  to  the 
convective  energy  transport  in  the  plume.  Nominally,  D  <  1.  As  the  fuel  feed 
velocity  gets  smaller  and  the  fire  base  size  gets  larger,  D  increases 
significantly  indicating  the  dominance  of  radiant  energy  loss.  The  parameter 
C,  evidently,  is  a  modified  Froude  number,  »  1.  Morton,  thus,  appears  to 
cover  the  important  problem  of  opaque  strongly  buoyant  plume  in  distinct 
contrast  with,  and  dramatic  complement  to,  the  transparent  weakly  buoyant 
plume  scaling  suggested  by  Murgai . 

9 . 3  Induced  Wind  Flow:  Some  Estimates  and  Some  Measurements:  An 

examination  of  the  scaling  rules  developed  in  section  (a)  above  indicates  that 
the  far-field  entrainment  velocity  ue  depends  on  height  z  above  the  source  of 
strength  v,  as 

1/3 

u  z 

—  »  constant,  near  0.66  (9.26) 

Q 

The  constant  itself  may  be  somewhat  smaller  in  tho  near-field  due  to  the 
observation  of  Morton  that  for  strongly  buoyant  plumes,  a=  ( p/  r>-»)  1  / 2  .  cw. 


(9.23) 
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Priestley's  Wmax  =  Wmax(wo«zo^z)  relation  can  be  transformed  to  the 
following  relation  for  entrainment  velocity. 

1/3 


V 


(  25  a\i/3  f 

/  Zn\21 
1  (  0  1 

w  )  L 

1  "l  r) 

(9.2 


Assuming  a  -  0.1,  the  plume  spreads  with  a  slope  of  0.1  so  that  for  a  source 
of  radius  b0,  the  virtual  source  lies  at  a  distance  of  zQ  =  10xbo  below  the 
real  source.  While  the  far-field  solution  Indicates  that  as  z  z0,  ue  tends 
to  a  large  value,  the  Priestley  solution  implies  that  ue(z  +  z0)  +  o  which  is 
physically  quite  reasonable. 

Taylor  [44]  discusses,  in  rather  general  terms,  the  Issue  of  air  Inflow 
into  fire  plumes  to  orimarily  highlight  the,  now  well-known,  fact  that  both 
the  axial  velocity  w  and  lateral  (entrainment)  velocity  ue  vary  with  height  z 
raised  to  a  power  of  -1/3.  This  recognition  of  the  proportionality  ue  «w, 
which  Is  a  latent  consequence  of  the  involved  turbulent  mixing  mechanism,  Is 
indeed  profound,  for  the  seed  is  sown  to  the  eventual  development  of  the 
concept  of  entrainment  constant. 

Faure  [50]  points  out  that  in  the  great  'Hamburg  raid,  two  out  of  every 
three  buildings  over  an  area  of  11  km2  were  ablaze  just  within  20  minutes 
after  the  first  wave  of  bombers  had  passed.  A  3  km  high,  2  1/2  km  base 
diameter  column  of  burning  gases  was  noted.  At  a  distance  of  2  1/2  km  from 
the  fire  zone,  winds  as  high  as  17  to  53  km/hr  were  induced  to  uproot  trees 
0.9  m  in  trunk  diameter.  Due  to  a  40  km?  fire  in  Tokyo  in  1945,  induced  winds 
were  noted  to  be  45-88  km/hr  at  a  distance  of  1  1/2  km  from  the  periphery  of 
the  fire;  upward  flow  was  so  strong  that  aircraft  flying  at  an  altitude  of  1 
800  meters  over  the  fire  zone  were  overturned. 

Faure  reports  an  important  experimental  fire  on  a  325-hectare  (803-acre) 
plot  of  1  m  high  bushes  and  thornbroom.  The  flames  were  observed  to  be  as 
high  as  10  to  12  m.  The  5  to  8  m/s  prevailing  wind  did  not  seem  to  have  much 
of  an  effect  on  the  smoke  column  upto  about  500-600  m  altitude.  The  fuel  load 
was  estimated  to  be  near  1.7  kg/m2  (1/3  lbm/ft2)  and  the  fuel  element 
combustion  time  was  about  35  minutes.  Assuming  that  the  fuel  has  a  heat  of 
combustion  of  about  16,000  kJ/kg,  the  heat  release  rate  is  estimated  to  be 
about  3.7  x  IQ®  megawatts.  Also  assuming  stoichiometric  air/fuel  ratio  to  be 
3.3  m^/kg  (53  ft3/lbm),  the  air  consumed  in  combustion  is  expected  to  be  1.6  x 
10®  m3  total  at  a  rate  of  8,000  m®/s.  The  average  required  entrainment 


velocity  over  a  rim  of  10  m  high  around  the  periphery  of  equivalent  diameter 
of  2,000  m  thus  is  about  1.25  m/s.  The  observed  air  speed  100  m  from  the  fire 
ranged  between  10  to  15  m/s  which  at  the  fire  rim  Is  expected  to  reflect  In  an 
actual  entrainment  speed  of  about  11-17  m/s.  Thus,  the  air  actually  entrained 
Is  over  10  times  larger  than  that  required  for  combustion;  this  Is  in  total 
agreement  with  the  general  knowledge  that  nearly  all  fires  Induce  a  1,000% 
excess  air.  With  this  Indicated  excess  air,  the  fire  gases  entering  the  base 
of  the  plume  can  be  estimated  to  be  at  about  400  K.  The  'smoke  column'  was 
noted  in  the  test  to  have  reached  an  altitude  of  about  2300  m  where  it  became 
horizontal . 

Faure  also  attempts  to  draw  the  scale  modeling  rules  for  the  flow  and 
combustion  phenomena  associated  with  very  large  fires  and  the  plumes  over 
them.  But  little  is  done  to  demonstrate  the  utility  of  the  developed  rules. 

A  similar  comment  is  valid  also  for  an  elegant  note  prepared  by  Byram  [51], 

An  investigation  of  the  turbulent  natural  convection  plume  above  a 
circular  ring  fire  has  been  reported  by  Lee  and  Ling  [52].  As  a  fire  grows 
gradually,  say  in  a  circular  pattern,  the  interior  of  this  circular  area  will 
be  'burnt  out',  sooner  or  later.  A  result  Is  a  fire  which  is  a  ring  of  finite 
width  spreading  outwards.  A  dimensional  analysis  is  presented  to  arrive  at 
the  following  velocity  and  temperature  scaling  in  terms  of  our  earlier-defined 
characteristics . 
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where  w  is  the  upward  velocity,  u  is  the  inward,  radial,  (entrainment) 
velocity,  r0  is  the  mean  radius  of  the  ring-fi re-base,  -96/71  Is  the  Rouse 
constant  and  Q  is  the  buoyancy  input  [r02  umax  g(pat-  p)max/  pj(m4/s3) . 

These  scaling  relations  are  essentially  same  as  those  of  Priestley  and  Ball 
and  of  Morton,  Taylor  and  Turner.  Measurements  made  in  the  plume  of  a  propane 
fire  over  a  29-inch  diameter  burner  have  been  correlated  on  the  basis  of  the 
developed  dimensional  analysis.  Among  a  number  of  interesting  conclusions, 
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the  following  two  are  of  our  immediate  consequence.  First,  as  one  plots  the 
entrainment  velocity  U  as  a  function  of  nondimensional  radius  (R  =  r/r0),  for 
various  axial  locations  Z(  =  z/r0)  and  various  characteristic  velocities 
u)(  =  [ 6^q/ ( 2 r0 (  +  l))]1/3,  one  finds  that  It:  (a)  is  virtually  independent 

of  u;  (b)  Increases  with  increasing  Z  (Umax  -  2  at  Z  «  0.5  and  «  2.5  at  Z  « 

4);  and  (c)  increases  from  zero  at  R  *  0  to  a  maximum  near  about  R  =  0.5  and 
then  decays  as  R  increases;  at  the  rim,  R  -  1,  U  is  about  50%  of  Its  maximum 
value.  In  physical  terms,  the  maximum  entrainment  velocity  noted  in  these 
29-inch  diameter  fires  is  about  6.5  ft/s.  That  the  maximum  entrainment 
velocity  occurs  within  the  fire  with  the  attendant  pressure-field  consequences 
and  that  the  velocity  and  temperature  fields  scale  well  within  the  context  of 
a)  are  perhaps  the  two  most  Important  contributions  of  this  work  of  Lee  and 
Ling. 

W.  J.  Parker,  Corlett  and  R.  T.  Lee  [53]  conducted  scale  model  tests  of  a 
mass  fire  (FI ambeau^ Countryman  [54])  using  a  0.61  m  x  0.61  in  (2  ft  x  2  ft) 
electrical  heat  source  of  strength  12.7  kW/m2.  The  Flambeau  fire  was  335.3  m 
x  320  m  (1100  ft  x  1050  ft)  in  extent.  The  length  ratio  between  the  Flambeau 
and  the  model  thus  is  537.  At  18  minutes  after  ignition,  the  Flambeau  test 
indicated  a  heat  release  rate  of  about  200  kW/m2  (1060  Btu/ft2  min).  At  an 
elevation  of  15.24  m  (50  ft)  at  the  center  of  one  edge  of  the  Flambeau  plot, 
Induced  velocity  was  measured  to  be  fluctuating  between  about  3  and  5  m/s  (10 
and  16  ft/s),  steady  after  about  18  minutes  past  ignition.  The  Parker, 
Corlett,  Lee  model  yields  a  velocity  of  0.195  m/s  at  a  height  of  0.028  m 

(which  is  homologous  to  50  ft).  With  the  scaling  argument  that  u  <*  A,  these 
modelers  convert  this  measured  model  velocity  to  the  prototype  Flambeau 

velocity  of  <iflal<bMU  -  Umode,  x  /(Ln  OTlbeau/l-11Iode1 )  -  0.115/OT  »  4.5?  m/s, 
which  is  within  the  range  of  the  Flambeau  fluctuation  measurement.  The 

scaling  argument  that  u  *  /L,  however,  is  in  some  variance  with  the  general 
knowledge  of  plume  scaling  discussed  earlier  in  this  chapter.  Some 
considerations  of  this  work  related  to  the  specific  details  of  the  combustion 
zone  will  be  discussed  later  in  section  (f)  of  this  report. 
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9.4  Firebrand  Transport:  The  convective  currents  produced  by  the  fire 

will  bring  about  an  enhanced  air  supply  to  the  burning  fuel;  the  feed-back 
heat  flux  is  thereby  enhanced  as  well.  Intense  heating  is  expected  to  produce 
mechanical  disruption  of  the  burning  fuel  elements  to  generate  firebrands  — 
namely  fragments  of  fuel  which  continue  burning  inflight.  The  blast-produced 
debris  fuel  piles  are  expected  to  be  copious  sources  of  firebrands.  Given  the 
initial  size  distribution  of  the  firebrands  along  with  upward  convective  plume 
velocity  and  the  prevailing  wind,  the  problem  at  hand  is  to  map  out  the 
firebrand  trajectories.  Continued  pyrolysis  (and  burning)  of  a  firebrand  in 
the  course  of  its  flight  is  itself  expected  to  significantly  influence  the 
trajectory.  The  topics  of  concern,  in  order,  are  the  production,  consumption 
(life-time)  and  transport  of  the  brands.  There  exists  little  quantitative  or 
qualitative  information  on  the  issue  of  production.  One  can  note,  however, 
that  both  the  number  and  size  of  brands  produced  may  be  expected  to  depend 
upon  the  rate  of  heating  of  the  fuel  (and,  hence,  on  the  stage  of  growth  and 
intensity  of  burning  of  the  fire)  and  upon  the  nature  of  the  fuel  itself. 

There  is  urgent  need  to  develop  these  dependencies  which  are  perhaps  amenable 
only  to  empirical  studies. 

Combustive  consumption  of  wood  particles  in  the  forced  convection 
conditions  of  flight  is  described  by  Tarifa,  del  Notario  and  Moreno  [55], 

Both  the  density  of  the  solid  (py,*)  and  the  size  of  the  particle  (r)  are  noted 
to  decrease  with  time  (t)  according  to 

fw  =  .  1  ,  (9.30) 
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where  the  subscript  o  represents  the  initial  conditions,  w  is  the  relative 
wind  speed  about  the  firebrand.  The  constants  n.  S  and  &  (with  units 
respectively  of  l/sec^,  m-/s  and  m)  are  empirically  determined  combustion 
properties  of  the  wood  with  a  given  moisture  content.  Tarifa  notes  that  the 
relative  wind  speed  w  is  not  significantly  different  from  the  particle's 
instantaneous  terminal  velocity  wf  which  depends  on  the  instantaneous  size  as 
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well  as  mass  of  the  particle  as 
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Thereby,  Tarlfa  obtains  the  following  implicit  relation  for  the  particle  size 
variation  with  time. 
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Wfo  in  these  relations  is  the  initial  terminal  velocity  of  the  brand  prior  to 
any  significant  reduction  in  size  and  mass  is  experienced. 

Describing  then  the  equation  of  motion  of  the  flying  brand  by  balancing 
its  Inertia  with  drag  and  gravity  forces,  Tarifa  computes  the  flight 
trajectory  numerically  under  conditions  of  constant  wind.  The  fall  velocity 
can  be  expressed  as 


wf(t)  = 


1  n»(t)  g  [I/2 
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where  m  is  instantaneous  mass  of  the  brand,  g  is  gravity,  Cq  is  drag 
coefficient,  Pa>  is  ambient  air  density  and  A  is  the  brand  frontal  area.  The 
main  difficulty,  even  if  the  brand  geometry  is  simplified  to  spheres  and  short 
cylinders,  was  that  both  the  drag  coefficient  and  the  mass  of  the  brand  were 
comollcated  functions  of  time  due  to  variation  in  size  and  density.  Among  the 
computed  quantities  is  the  time  (tg)  required  for  the  brand  to  reach  the 
ground.  Comparing  this  time  with  the  experimentally  determined  burn-out  time 
(tt>)  of  the  brand,  one  can  reach  certain  critical  conclusions.  If  tb  >  tg, 
then  firebrand  reaches  the  ground  while  still  burning,  thus  capable  of 
initiating  a  fire  upstream.  If  tb  <  tg,  the  brand  is  burnt-out  in  the  flight. 
The  maximum  horizontal  range  within  which  spotting  fires  can  be  initiated  by 
the  flying  brands  is  given  by  the  condition  tb  =  tg. 

Notwithstanding  the  need  to  employ  empiricism  in  describing  the 
combustion  of  the  firebrands,  it  is  fitting  to  note  that  the  concept 
underlying  the  firebrand  flight  methodology  of  Ref.  [55]  is  indeed 
meritorious.  Tarifa  himself  promplty  recognizes  this  and  presents  a 


dimensional  analysis  in  Ref.  [56],  Two  cases  are  considered:  first,  the 
firebrand  undergoes  combustion  in  a  wind  of  constant  speed;  and  second,  the 
firebrand  burns  as  it  flies  at  Its  final  fall  velocity. 

Denote  the  ratio  of  final  fall  velocity  of  the  brand  to  the  initial  fall 
velocity  wf/wr*0  by  the  symbol  5.  For  spherical  and  cylindrical  brands,  the 
experimental  data  indicate  that  this  ratio  depends  only  on  a  nondimensional 
quantity  Ui  defined  as 
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where  x  =  wt/d0  (with  wind  speed  w,  time  t  and  the  particle  diameter  d0); 

Re  =  ud0/v„  is  the  Reynolds  number  ( \>M  being  the  kinematic  viscosity  of 

ambient  air);  5  =  is  the  ratio  of  the  ambient  air  density  to  the 

initial  density  of  the  brand-wood;  and  X,  is  a  geometric  ratio  of  the  brand 

( Xl  =  length/diameter  ratio  for  cylinders  and  xi  3  '1  for  spheres).  The  burn 
out  time  X5  is  found  from  the  experimental  correlations  to  be  given  by 

Ul  =  Uib  =  2.6. 

In  the  second  case,  the  firebrands  are  considered  to  be  burning  within  an 
air  current  whose  velocity  is  constant  and  equal  to  the  brand's  final  fall 
velocity.  In  this  case  x  and  are  defined  on  the  basis  of  the  initial  fall 
velocity  wf0.  Then  the  expression  for  1)  remains  same  as  in  the  first  case  for 
cylindrical  and  spherical  particles.  For  square  plate  particles,  however, 

U  =  U]_  x  0.6x2  where  xj  represents  the  additional  required  length  ratio  in  the 
particle  description. 

In  order  to  further  employ  Tarifa's  methodology,  considerable  further 
work  on  the  behavior  of  firebrand  burning  is  required.  Various  particle 
shapes  and  ranges  of  initial  size  have  to  be  considered  on  different 
materials.  Two  recant  works  in  these  directions  are:  Lee  and  Heilman  [571 
and  Muraszew,  Fedele  and  Kuby  [58], 

9.5  Fire  and  Wind/Fire  Storms:  When  the  fire  is  large  enough  in  size 

and  in  intensity,  the  convective  flows  prompted  by  it  may  be  strong  enough  to 
noticeably  interfere  with  the  atmosphere  itself.  Considerable  research  in 
this  area  has  been  reported  in  the  literature  by  Morton  [59],  Long  [60], 
Muraszew,  Fedele  and  Kuby  [58],  Small  and  Erode  [61]  and  Carrier,  Fendell  and 
Feldman  [52], 


Morton's  paper  is  essentially  a  qualitative  description  of  the  flows 
Induced  by  large  fires,  conflagratory  fire  spread  as  influenced  by  the  ambient 
winds,  and  the  origin  of  blow-up  -(or  back-)  fires  which  generate  vorticity  to 
cause  fi re-tornadoes.  Two  points  made  in  this  paper  are  noteworthy. 

First,  whereas  the  entrainment  occurs  at  all  heights  with  a  velocity  about  an 
order  of  magnitude  smaller  than  the  mean  upward  velocity  in  the  plume,  the 
radial  pressure  field  near  the  ground  tends  to  amplify  the  local  entrainment. 
This  amplified,  near-ground,  entrainment  reflects  in  the  induced  winds. 
Furthermore,  it  is  expected  to  intensify  the  combustion  rate  of  the  loosely 
packed  fuel  bed  not  only  by  the  increased  air  supply  and  its  mixing  with  the 
fuel  vapor  but  also  by  augmenting  the  heat  flux  from  the  flame  to  the  fuel 
elements. 

The  second  point  concerns  with  the  effect  of  (back-)  fires  which  are 
started  by  such  agents  as  firebrands.  Known  also  as  blow-up  fires,  they  tend 
to  develop  strong  rotational  flow  with  a  concentrated  vortex  cone  to  yield  the 
phenomena  of  fire-whirls  and  fire-tornadoes.  These  phenomena  themselves 
serve  to  transport  the  firebrands  to  greater  altitudes  and  greater  horizontal 
reaches.  Thus,  the  problem  obviously  is  of  strong  nonlinearily  within  which 
discernment  of  cause  and  effect  is  not  at  all  a  clear-cut  exercise.  The 
strong  buoyant  flow  in  the  plume  and  the  associate entrainment  are  also  knowo 
to  amplify  the  vorticity  which  is  dilutely  present  in  the  atmosphere. 

Long  [60]  presents  an  analysis  of  the  fire-whirl.  The  buoyant  fire  plume 
entrains  and  amplifies  the  atmospheric  vorticity.  While  the  vorticity  due  to 
earth 1 s  rotation  is  of  the  order  of  10"4/sec,  a  hundred  times  larger  vorticity 
is  required  to  produce  a  fire-whirl.  Sources  of  such  a  larger  ambient 
vorticity  may  be  local  meteorological  eddies  or  cyclones  with  length  scales 
less  than  a  mile  or  so.  Gravitational  instability,  light  wind  and  low 
humidity  of  air  are  noted  by  Long  to  be  conditions  conducive  for  a  firestorm. 

Long  deduces  an  interesting  scaling  scheme  for  firestorms  (ignoring 
radiation,  treating  the  fire  as  a  buoyancy  source  by  simply  specifying  a 
temperature  Tf  in  the  fire  of  diameter  d,  higher  than  in  the  ambient  air  T., 
neglecting  the  Coriolis  force  due  to  earth's  rotation  (i.e.,  the  Rossby  number 
is  much  greater  than  unity),  considering  that  the  curvature  of  the  earth  is 
inconsequential  and  that  all  molecular  transport  processes  are  unimportant). 
The  nondimensional  parameters  to  be  preserved  between  the  model  and  prototype 
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o  =  adiabatic  exponent,  Cp/Cv  *<  1.4 
The  parameters  s,  m  and  a  are  thermodynamic  in  nature,  p  is  potential 
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density,  a  conserved  quantity,  of  the  fire  gases;  pQ  is  potential  density  of 
the  ambient  air;  a p  =  pQ  -  p.  If  the  model  and  prototype  fires  are  of  same 
temperature,  s  is  automatically  preserved.  So  also  is  a.  The  parameter  m  is 
indicative  of  the  compressibility  effects,  (g  is  gravity,  a  is  angular 


velocity  of  solidus  rotation  of  the  ambient  fluid;  d  is  fire  diameter;  T5*  is 
potential  pressure  of  the  atmosphere.)  Even  for  a  fire  as  large  as 
d  m  1500  m,  the  parameter  m2  is  seldom  larger  than  about  10“3  for  typical 


values  of  n,  and  Ap/p0*  It  could  thus  be  safely  disconcerning. 

Of  the  remaining  three  parameters,  e  is  a  measure  of  the  gradient  of 
potential  density  with  altitude.  In  its  definition,  ais  a  constant  related 
to  the  stability  of  the  atmosphere.  For  a  near-neutrally  stable  atmosphere 
which  is  recognized  to  be  conducive  for  a  firestorm,  o  «  o  so  that  »  o. 

Thus,  Long  rules  out  the  importance  of  e  in  the  scaling  endeavors  so  long  as 
care  is  taken  to  keep  it  near  zero.  One  then  is  left  with  only  two  parameters 
of  concern.  k‘  is  the  ratio  of  surface  roughness  length  k  (due  to  buildings 
in  a  city  or  trees  in  a  forest)  to  the  fire  diameter  d.  Keeping  k'  invariant, 
therefore,  implies  satisfying  geometric  similarity. 

Since  g  and  s  are  same  in  both  the  model  and  the  prototype,  one  has  to 
adjust  the  external  fluid  rotation  n  with  fire  size  d  so  as  to  maintain  the 
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important  parameter  F  Invariant.  Thus,  to  simulate  ad®  1500  m  fire  (in  an 
atmosphere  containing  an-  10-3/s  rotation)  by  a  model  of  d  s  15  m,  n  has  to 
be  adjusted  to  10-2/s.  This  can  be  achieved  by  mechanical  means. 

While  k'  and  F  thus  are  the  primary  scale  parameters  to  simulate  a 
firestorm,  Long's  analysis  indicates  that  the  pressure,  density,  radial 
location,  velocities  and  time  would  scale  according  to  the  following 

nondimensionalization  with  L  =  [n^dVgf Ap/"p0) ] 1  3. 

P*  =  "P"’ / gL Ap  ;  P*  a  p'/a7  ;  r*  =  r/L  ;  (9.42 

v*  =  v/[gLA7/70l1/2  and  t*  a  (t/DCgLi^l 1/2  .  (9.4: 

If  there  were  no  rotation  whatsoever  Is  present,  n  will  be  zero;  the 
characteristic  length  becomes  d  in  these  definitions;  F  will  be  zero,  m2  and 

e  respectively  become  gd/CF  /"p  )  and  ad/(Ap/p  );  m  will  still  be  negligible; 
e  will  remain  near  zero  for  both  the  model  and  prototype;  and  the  modeling 
endeavors  reduce  to  keeping  thermodynamics  invariant  and  preserving  geometric 
similarity.  Suffice  it  to  note  in  conclusion  that  Long's  proer  holds  a  wealth 
of  ideas  pertaining  to  scaling  mass  fires,  especially  firesorms;  it  should 
be  a  standard  reference  in  all  fire-modeling  efforts. 

9.6  The  Combustion  Zone  Specifics:  Once  the  heat  output  of  a  fire  over 

a  given  area  at  altitude  zero  is  given,  it  is  clear  from  the  preceeding 
discussions  that  the  character  of  buoyant  plume  flow  is  scaleahle.  To 
determine  the  pyrolysis  and  burning  rate  of  the  fuel  deris  is  thus  a  crucial 
aspect  to  be  discussed  below. 

As  a  condensed-phase  fuel  element  is  heated,  convectively  due  to  contact 
witn  the  hot  flame  gases  or  radiatively  due  to  remote  flames  and  glowing 
charcoal,  It  first  heats-up  and  then  undergoes  pyrolysis  (or  simple  physical 
vaporization)  to  yield  combustible  pyrolyzate  gases  and  vapors.  Upon  mixing 
with  the  oxygen  of  Hie  ambient  air,  these  pyrolyzates  result  in  fresh  flames 
which  may  themselves  furnish  further  fuel-heating  energy.  The  problem  of 
combustion  intensity  thus  entails  a  study  of  heat  and  mass  transfer  across  the 
fuel  element  boundary.  While  conductive  heating  and  subsurface  pyrolysis 
processes  predominate  within  the  fuel  element,  convection  and  oxidation 
reactions  dominate  in  the  gas-phase  (and  at  the  interface).  Thermal  radiative 
heat  exchange  may  become  relevant  both  within  and  outside  the  porous  fuel 
element. 


When  an  Isolated  fuel  element  Is  subjected  to  radiative  and  convective 
heating,  the  steady  state  pyrolyzate  efflux  rate  m"  (kg/m^s)  can  he  written, 
by  a  linear  approximation,  as 

*"  ■  «'r,d  +  h<T.  '  V^sg  ,9'44) 

where  q'Vad  is  radiant  flux  received  by  the  surface,  h  is  heat  transfer 
coefficient,  T.  is  the  flame  gas  temperature,  Tw  Is  the  surface  temperature  of 
the  pyrolyzing  fuel  element  and  Lsg  is  the  latent  heat  of  pyrolysis  (J/kg 
pyrolyzate).  With  Cpg,  Kg  and  i  respectively  representing  the  gas-phase 
specific  heat,  conductivity  and  the  specimen  characteristic  dimension,  the 
above  equation  can  be  written  nondimensi onally  as  follows: 
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The  left-hand-term  is  nonaimensional  pyrolyzate  efflux  rate.  The  first 
term  on  RHS  accounts  for  the  radiant  heating  while  the  second  term  accounts 
for  convective  heating.  The  quantity  (hz/Kg)  is  known  as  Nusselt  number,  a 
function  of  fluid  mechanics  of  flow  (i.e.,  of  Reynolds  or  Grashof  number  and 
Prandtl  number).  The  enthalpy  ratio  [CpgCT*  -  Tw)/LSg]  is  a  descriptor  of  the 
thermodynamics  of  combustion  (numerator)  and  pyrolysis  (denominator);  this 
ratio  is  known  as  the  B-number. 

The  role  of  radiative  heating  may  be  viewed  in  different  view-points. 
First,  rewriting  the  preceeding  relation  in  the  following  form. 
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the  effect  of  radiation  may  be  taken  as  an  apparent  increase  in  the  flame 
temperature  to  augment  the  convective  heat  transfer  rate.  Alternatively, 


(9.47) 
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whereby,  rad i 1 1  on  may  be  viewed  to  have  increased  the  convective  heat  transfer 
coefficient.  Finally,  if  q"rad  /[h(Tw  *  Tw)]  is  small  compared  to  unity, 


f”"CP9‘' 

iu' 

r  (T  -  T  ) 
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K9  J 
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L  (1  -  [q“  ,  ./h(T  -  T  )]) 

L  sg  ^  rad  »  w  „ 

Implying  that  the  presence  of  radiant  heating  tends  to  manifest  in  a  reduction 
of  the  enthalpy  of  pyrolysis. 

This  simple  linear  analysis  is  thus  quite  instructional  as  to  the  manner 
in  which  thermal  radiation  Intrudes  to  augment  the  convective  pyrolysis  rate 
of  a  fuel  element.  Exact  analyses  solving  the  differential  equations 
governing  the  heat  and  mass  transfer  process  are  available  to  remove  the 
present  simplification  of  linearity.  But  the  fundamental  nondimensional 
parameters  remain  unaltered. 

It  then  remains,  in  scaling  the  heating  process  of  a  fuel  element,  a 

trivial  task  to  estimate  the  radiant  flux  q"rad  and  to  obtain  the  Nusselt 
number  (h Jt/Kg)  as  a  function  of  Reynolds  (or  Grashof)  and  Prandtl  numbers. 
Existing  heat  transfer  literature  enables  us  to  perform  this  task  quite  well. 
The  main  difficulty,  however,  lies  in  the  fact  that  the  fires  of  present 
concern  involve  fuel  beds  which  are  nonuniform  and  fuel  elements  in  these  beds 
are  random  in  size,  shape  and  orientation.  This  difficulty  may  be  overcome, 
perhaps,  only  by  a  stochastic  description  of  the  structure  of  the  fuel  bed. 
Such  descriptions,  however,  are  presently  unavailable.  Not  withstanding  this 
unavailability,  it  is  perhaps  worth  noticing  how  Rothermel  coped  with  this 
fuel  description  problem. 

In  order  to  study  the  forced  convective  combustion  of  fuel  element  debris 
piles,  Strasser  and  Grumer  [63]  reported  an  interesting  experiment.  4  kg  of 
2x2x4  cm  wood  pieces  were  piled  randomly  in  a  wire  basket,  ignited  with  a 
methanol  accelerant  and  combustion  rate  was  monitored.  The  rate  quickly 
attains  a  maximum  (»  4.3  g/s)  and  then  gradually  decays.  The  air  flow  induced 
into  the  fuel-bed  was  modeled  by  a  well-known  relation  of  the  pressure  drop 
through  packed  beds  with  turbulent  flow  [64],  In  the  balance  of  the  pressure 
drop  with  buoyancy  and  inertia,  it  has  been  found  that  the  inertia  effect  is 
minimal.  The  velocity  of  air  induced  through  the  peripheral  area  of  the 


vertical  circular  cylinder  of  a  debris  pile  (height  a  Ht,ed»  diameter  a  d^ed) 
was  found  to  depend  upon  the  ambient  air  temperature  T*,  the  burned  gas 
temperature  exiting  out  of  the  pile  at  the  top,  the  diameter  of  the  fuel 
particles  dpart-)c1e  composing  the  bed  whose  void  fraction  is  5,  the  friction 
factor  f  which  Is  a  weak  function  of  the  Reynolds  number,  and  a  particle  shape 


factor  x. 
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(5.49) 


The  velocity  thus  calculated  was  not  only  in  excellent  agreement  with  the 
measured  but  also  equal  to  the  stoichiometric  requirement  for  the 
wood-burning.  Recall,  from  our  estimates  In  the  context  of  Faure's  field 
experiment,  that  wood-like  fuels  require  about  3.3  m^  («4  kg)  of  ambient  air 

to  stolchlometrlcal ly  burn  1  kg  of  the  fuel.  Thus,  the  burning  rate  try  (kg 
fuel/sec)  of  the  debris  basket  appears  to  obey  the  relation 
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[9.50) 


While,  thus,  the  air  Induced  into  the  fuel  bed  Itself  satisfies  the 
stoichiometry,  the  flames  are  expected  to  project  above  the  fuel  bed  to  a 
height  of  the  order  of  fuel  bed  diameter  for  small  but  turbulent  fires  and  of 
the  fuel  bed  height  for  large  fires.  It  Is  these  flames  which  induce  great 
quantities  of  excess  air  to  dilute  the  flame  gases. 

The  maximum  burning  rate  m_,„  measured  by  Strasser  and  Grumpr  is  such 

iT)d  X 

that  the  burn  time  of  a  fuel  element  of  initial  mass  of  m0  is  in  the  range 

1200  <  ^o/%ax  <  2800  sec  for  Douglas  fir  sticks  In  the  range  2  cm  <  d  part^cie 
<  9  cm.  When  a  fan  blows  air  into  the  flame  column,  no  discernible  change  In 
burning  rate  was  noticed.  However,  if  it  blows  air  Into  the  fuel  bed,  the 
burning  rate  is  enhanced  somewhat  in  proportion  to  /( uwl nd/dpa rt i cl e ^  * 

Attempts  to  progress! vely  increase  the  size  of  laboratory  fires  have 
always  yielded  the  result  of  disintegration  of  the  flame  into  multiple 
flamelets.  (See  also  Corlett  [65],  Blackshear  [66]  and  Wood,  Blackshear  and 


Eckert  [67].)  When  the  fire  base  diameter  Is  less  than  10  m,  the  visible  flame 
height  Is  approximately  equal  to  the  base  diameter.  At  larger  base  diameters, 
however,  the  flame  height  Is  much  smaller  than  the  base  diameter  and  the  flame 
appeared  to  be  composed  of  multiple,  Individual  flamelets  which  move  around  the 
base  in  an  unsteady  way  almost  as  though  the  heat  feed-back  and  vaporization 
processes  are  nonuniform. 

Attempts  to  model  the  Interaction  of  multiple  adjacent  flames  (Putnam  and 
Spelch  [68],  Welker  and  Sliepcevlch  [31])  are  less  than  completely  successful. 
Whereas  flames  separted  widely  tended  to  mutually  augment  the  burning  rate 
(presumably  due  to  enhanced  heat  feed-back),  flames  separated  but  only  narrowly 
tended  to  experience  lowered  burning  rate  (presumably  due  to  competition  for 
oxygen).  Reference  [63]  presents  some  data  relevant  to  this  behavior.  When 
two  26  cm  high  fuel  beds  composed  of  4  cm  cubes  of  Douglas  fir  fuel  elements 
were  burned  side-by-side,  the  maximum  burning  rate  Increased  from  3.3  g/s  when 
separation  Is  large  to  a  maximum  of  about  5  g/s  at  a  separation  of  about  0.1  m 
and  then  decreased  to  about  3.7  g/s  when  separation  is  zero.  It  was  concluded 
that  radiative  heat  transfer  between  proximate  fires  augment  the  burning  rate 
less  than  convective  heat  transfer  does.  Strasser  and  Grumer  [63]  arrive  al: 
the  opinion  that  a  metastable  atmosphere  Is  perhaps  a  prerequisite  for  multiple 
fires  to  merge  to  form  a  firestorm  and  that  in  such  a  situation,  combustion 
occurs  within  (and/or  very  close  to)  the  fuel  bed. 

Parker,  Corlett  and  B,  T.  Lee  [53]  also  come  to  a  similar  conclusion  based 
on  the  following  two  points.  First,  in  large  fires,  molecular  transport 
processes  are  almost  insignificant  relative  to  turbulent  mixing.  The  Grashof 
number  is  thus  realized  to  be  so  large  as  to  be  irrelevant  In  mass  fire 
scaling.  Second,  the  flames  and  product  gases  in  large  fires  are  optically 
thick  in  the  wavelength  range  of  concern;  the  Rosseland  diffusivity  of  such  a 
medium  is  infinitecimally  small  compared  to  the  turbulent  eddy  diffusivity  for 
heat  transport.  The  net  consequence  of  this  is  that  radiative  transport  within 
a  mass  fire  is  expected  to  be  dominant  only  in  a  thin  zone  near  the  fuel  bed. 

It  must  thus  be  permissible  to  model  mass  fires  by  employing  an  electrical 
(area)  heat  source  near  the  ground  surface  in  the  laboratory. 

Having  thus  simplified,  the  mass  fire  problem  is  reduced  to  a  thermally 
driven  fluid  mechanics  problem  whose  working  fluid  is  normal  air.  The 
combustive  energy  input  then  can  be  expressed  as  a  merely  nondimensional  heat 
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where  (g  i) 1/2  is  the  characteristic  velocity,  i  is  the  fuel  bed  dimension  and 

q"  is  the  heat  flux  (W/m^).  The  heat  release  rate  thus  being  proportional  to 

/I,  a  600  m  (■  2000  ft)  fire  releasing  energy  at  a  rate  of  about  3  MW/m^  ( «.  10b 
Btu/ft2/hr)  is  expected  to  be  scaleable  by  a  0.6  m  »  2  ft)  laboratory  fire 
with  a  heat  flux  of  95  kW/m2  («  30,000  8tu/ft2/hr).  This  small  a  laboratory 
fire,  albeit  turbulent,  is  indeed  a  weak  one  exhibiting  an  irregular  flame 
hovering  erratically  about  the  fuel  surface  ‘with  much  openness'  [53]. 

Corlett,  too,  concludes  that  without  the  effects  of  vertical  atmospheric 
vortlcity,  a  firestorm  will  most  probably  not  arise. 

While  these  considerations  attempt  to  capture  the  essential  heat  transfer 
mechanisms  between  the  flames  (and  hot  gases)  and  the  surface  of  a  fuel 
element,  it  is  fitting  to  briefly  examine  the  scaling  principles  involved  in 
the  manner  in  which  the  solid  fuel  element  utilizes  the  energy  received.  As 
the  solid  heats  up  due  to  transient  conduction,  it  eventually  undergoes  thermal 
degradation  at  all  depths.  The  pyrolyzates  thus  produced  transpire  out  of  the 
char  into  the  flaming  surroundings. 

Tarifa  [55]  notes  from  wind-tunnel  experiments  that  the  density  of  wood 
decreases  with  time  according  to  an  inverse  quadratic  function  of  time: 

0w/pwo  =  1 / [1  +  n  t2]  where  n  is  an  empirical  constant  characteristic  of  the 
nature  of  wood.  He  also  finds  that  the  diameter  of  the  fire  brand  decreases 
according  to  a  square-law:  (d/d0)2  =  1  -  [4(g  +  ws)/d02]t  which  is  analogous 
to  the  well-known  square-law  associated  with  liquid  fuel  droplet  combustion. 
(Recall  that  3  and  6  are  empirical  constants  and  w  is  wind  speed.) 

An  analysis  of  charring  conduction  presented  by  the  present  author  [69] 
sheds  further  light  on  the  nature  of  combustion  of  wood-like  fuel  elements. 

The  main  result  of  this  analysis  is  the  following  time-dependency  of  the  mass  m 
of  the  charring  element. 


sill  =  *£.  »  /,P°  I..pc\  A  ■  V\(Jn)/2  (9.52) 
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where  m0  Is  Initial  mass,  pc  and  p0  are  final  char  and  initial  wood  densities, 
t  is  time  and  the  Index  j  is  o  for  an  infinite  slab,  1  for  a  long  cylinder  and 
2  for  a  sphere.  The  characteristic  burn-out  time  t*  is  given  by 
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where  i0  is  the  initial  half-thickness  of  the  particle,  LSg  is  enthalpy  of 
pyrolysis  (kJ/kg  pyrolyzate),  Cp0,  cxo  are  respectively  the  specific  heat  and 
thermal  diffusivity  of  wood,  Tp‘  is  a  temperature  characteristic  of  the 

pyrolysis  process,  T0  is  wood's  Initial  temperature,  and  q"  is  the  heat  flux 
received  by  the  surface  from  the  ambient  gas-phase.  Two  major  points  are 
evident  from  these  results:  (a)  The  mass  loss  rate  depends  on  the  time 
function  (1  -t/t*)  raised  to  a  power  of  -1/2  for  a  slab,  0  for  a  cylinder  and 
+1/2  for  a  sphere.  For  cylindrical  fuel  elements,  thus,  the  rate  is  apparently 
independent  of  time,  (b)  The  terms  on  the  RHS  of  the  (c^t*/^2)  equation 
respectively  represent  the  thermochemistry  of  pyrolysis  and  thermal  diffusion. 
When  thin  fuel  elements  are  considered,  the  thermochemical  term  dominates  to 
make  visible  the  influence  of  pyrolysis  endothermicity  Lsg,  pyrolysis  kinetics 

Tp‘  and  the  heating  rate  q".  When  the  fuel  elements  are  thick,  however,  the 
second  term  dominates  indicating  that  thermal  diffusion  overwhelms  all  the 
thermochemical  and  thermokinetic  details  as  well  the  exposure  flux 
specification  itself.  With  lSg  «  300  kJ/kg,  pq  «  600  kg/m3,  K  »  150  kg/m3 
and  oq  »  10“7  m2/s,  the  above  result  becomes 
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where  t0  is  in  meters,  q"  is  in  Watts/m2  and  t*  is  in  seconds.  If  q"  is 
20,000  W/m2,  for  cylindrical  dowels,  t*  »  6.75  x  10^  ^  +  1.25  x  10® 
so  that  fuel  elements  thicker  than  about  a  centimeter,  or  so,  will  behave  as 
though  they  are  thermally  thick  and  the  pyrolytic  time  is  determined  more  by 
the  diffusion  effects  than  by  the  chemical,  kinetic,  and  exposure  effects. 
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SECTION  10 


CONCLUDING  REMARKS 

Based  on  the  foregoing  discussions,  It  is  eminently  clear  that  much  can  be 
done  to  Improve  our  understanding  of  the  behavior  of  a  large  real  fire  by 
employing  the  technique  of  scale  modeling.  At  the  heart  of  the  design  of  a 
scale  model  lies  the  need  to  recognize  the  physicochemical  mechanisms  involved 
in  the  character  of  the  modeled  fire.  This  recognition  can  evolve  from  either 
a  sophisticated  theoretical  formulation  of  the  problem  (but  with  no  need  to 
solve  it)  or  a  mere  dimensional  analysis  (in  which  a  dimensional  assessment  is 
made  of  all  the  Involved  dependent  and  independent  variables).  The  main 
objective  of  a  scale  model  design  exercise  is  to  identify  combinations  of 
variables  of  the  problem  so  that  a  relationship  among  these  combinations  may  be 
easier  determined  experimentally  without  destroying  the  generality  of  the 
relationship.  In  other  words,  the  behavior  patterns  of  the  phenomenon  under 
study  will  be  quantitatively  expressible  as  relations  between  the  variable 
combinations  such  that  the  proportionality  constants  and  indices  of  powers  are 
Independent  of  the  system  and  even  of  the  materials  employed.  The  combinations 
of  variables  are  nondimensional ;  they  serve  to  relate  the  scale  model 
homologously  with  the  prototype  and,  as  such,  are  known  as  scaling 
prescriptions. 

Becoming  thus  impressed  with  the  tremendous  proven  power  of  scale  modeling 
in  various  engineering  disciplines,  attempts  to  apply  this  art  in  studies  of 
fire  behavior  quickly  lead  one  to  conclude  that  complete  modeling  of  such  a 
complex  system  is  Impossible,  for  the  number  of  prescription  rules  far  exceed 
the  experimenter's  degrees  of  freedom.  The  reader  is  encouraged  to  refer  to 
Professor  Spalding's  superb  expose  [8]  to  truely  appreciate  the  difficulties 
associated  with  complete  modeling  of  combustion  systems. 

We  are  thus  convinced  that  we  should  become  content  with  modeling  of  the 
mass  fire  in  parts.  The  contents  of  the  present  report  indicate  that  the 
mechanistic  description  of  various  aspects  of  a  mass  fire  is  not  only  possible 
but  also  mostly  available  in  the  existing  literature.  These  aspects  range  from 
the  inception  of  the  fire  to  the  fully  developed  plumes.  The  state-of-the-art 
of  current  understanding  of  these  aspects  is  so  advanced  that  the  scaling  rules 
can  be  relatively  easily  deduced. 

The  question  facing  immediately  is:  Which  aspect  or  aspects  of  the  mass 


fire  Is  It  that  we  wish  to  scale  model?  The  answer  obviously  depends  upon  the 
"interests  of  the  experimenter  and  the  accuracy  and  urgency  of  the  required 
prediction"  [8].  Once  this  choice  for  partial  modeling  is  made,  scaling 
prescriptions  can  be  deduced  to  specifically  highlight  the  effects  sought; 
scale  model  experiments  can  be  designed  to  satisfy  these  prescriptions; 
experiments  can  be  conducted;  and  their  results  can  be  generalized. 

Evident  from  the  work  reported  here  are  some  areas  in  which  the  existing 
mechanistic  understanding  leaves  further  progress  desired.  Some  of  these 
areas  are:  the  mechanisms  of  fire  growth  within  enclosures  to  flashover;  the 
role  played  by  thermal  radiation  heat  transfer  In  a  number  of  contexts  of  fire 
growth,  conflagratory  speed,  intensity  of  burning  and  others;  the  mechanism  of 
production  of  firebrands  in  large  fires;  the  phenomenon  of  firebrand  burning 
Inflight;  and  the  specific  details  of  combustion  of  fuel  elements  in  the 
Intense  environment  of  debris  piles. 
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APPENDIX 


GLOSSARY  OF  TERMS 

The  multitude  of  fire  studies,  from  which  scaling  relationships  are  drawn 
here,  frustrate  all  deliberate  attempts  to  unify  the  nomenclature.  A 
section-by-section  nomenclature  appears  to  be  desirable  to  avoid  ambiguities 
notwithstanding  some  certain  redundancy. 


Section  5: 

c  a  constant  coefficient,  Table  I. 

D  distance  from  ground-zero,  m, 

Etotal  thermal  radiation  yield,  J. 
f  thermal  radiant  energy  partition  =  Etotal • 
F  a  function. 

A 

F  fluence,  J/m2. 

H  height  of  burst,  m. 

I  irradiance  at  the  target,  W/m2. 

m  an  index.  Table  I. 

n  an  index.  Table  I. 

P  thermal  energy  release  rate,  W. 

R  radius,  m. 

t  t i me ,  s . 

T  temperature  of  the  fireball,  K. 

W  weapon  yield,  J. 

<  absorption  coefficient  of  air,  1/m. 

p  density  of  rir,  kg/m3. 

p0  density  of  air  at  sea  level,  kg/m3, 
o  Stefan-Bol tzmann  Constant,  W/m2<4> 

t  transmissivity  of  air  mass. 
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Section  6: 


specific  heat  of  the  solid,  J./kg  K. 
gravitational  acceleration,  m/s2, 
convective  heat  transfer  coefficient,  W/m2K. 
height  of  the  target  specimen,  m. 
irradiance  on  target  surface,  W/m2. 
thermal  conductivity  of  the  solid,  W/mK. 
solid  target  specimen  thickness,  m. 
time,  s. 

Initial  temperature,  K. 
pyrolysis  temperature,  K. 
ambient  gas  temperature,  K. 

Ignition  temperature,  K. 
solid  thermal  diffusivity,  m2/s. 
solid  surface  emissivity. 
gas-phase  kinematic  viscosity,  m2/s. 
solid  density,  kg/m2. 


Section  7: 


gas-phase  specific  heat,  d/kg  K. 

activation  energy  of  oxidation  reaction,  J/kgmole. 

enthalpy  of  combustion,  0/kg. 

preexponential  factor  of  oxidation  reaction,  1/s. 

fuel  bed  length  along  the  flow  direction,  m. 

Oamkohler  number. 

energetic  potency  parameter. 

universal  gas  constant,  J/kgmole  K. 
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u  blast  wind  peak  velocity,  m/s. 

Yp -(  fuel  vapor  mass  fraction  at  the  fuel  surface. 

Section  8 
Subsection  3,1: 

B-j  Biot  number. 

CpS  solid  specific  heat,  J/kg  K. 

0  diffusivity  of  gases  through  porous  solid,  m2/s 

1l  nondimensional  diffusivity  0. 

E  activation  energy  of  smolder  reaction,  J/kgmole 

f  stoichiometric  fuel/oxygen  mass  ratio,  kg/kg. 

h  convective  heat  loss  coefficient,  W/m2K. 

hc  enthalpy  of  smoldering  combustion,  .)/kgmole. 

K,KS  solid  thermal  conductivity,  W/mK. 

Ji.is  solid  half-thickness,  m. 

R  universal  gas  constant,  J/kgmoleK. 

t  time,  s. 

Tp  pyrolysis  temperature,  K, 

u  smolder  spread  velocity,  m/s. 

V  nondimensional  u. 

Y0„  oxygen  mass  fraction  in  ambient  gas. 

Z  characteristic  smolder  reaction  rate,  kgmoie/m3 

5  a  Damkohler  number. 

pc  char  density,  kg/m3. 

pg  gas  density,  kg/m3. 

ps  virgin  solid  density,  kg/m^. 

q*  nondimensional  ambient  temperature. 
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Subsection  8.2: 


cpg 

gas  phase  specific  heat,  <3/kg  K. 

cps 

solid  specific  heat,  0/kg  K. 

h 

gas  phase  convective  heat  transfer  coefficient, 

W/m2K. 

Ah  5 

solid  enthalpy  rise  from  Initial  to  Ignition  state,  *3 /kg 

Kg 

gas  thermal  conductivity,  W/mK. 

KSy 

solid  conductivity  in  normal  direction,  W/mK. 

l9 

gas  phase  conduction  length,  m. 

Is 

thin  fuel  bed  thickness,  m. 

n 

characteristic  length  of  relevance  to  Ri,  m. 

q" 

heat  flux  across  the  surface  of  flame  inception 

,  W/m2. 

q" 

g 

gas  phase  conduction  flux,  W/m2. 

Rl,R2 

upstream  and  downstream  radiant  fluxes,  W/m2. 

t 

time,  s, 

Tf 

flame  temperature,  K. 

Tvap 

vaporization  or  pyrolysis  temperature,  K. 

ambient  gas  temperature,  K. 

u 

flame  spread  velocity,  m/s. 

Ug 

opposing  ambient  wind  speed,  m/s. 

2 

a  dummy  parameter  of  relevance  to  Rj. 

OS 

solid  thermal  diffusivity,  m2/s. 

pg 

ambient  gas  density,  kg/m2. 

os 

fuel  bed  solid  density,  kg/m2. 

Subsection 

8.3: 

8 


transfer  number. 


gas  specific  heat,  J/kg  K. 
solid  specific  heat,  J/kg  K. 
stoichiometric  fuel/oxygen  mass  ratio,  kg/kg. 
gravitational  force  per  unit  mass,  m/s^. 

Galileo  number. 

heat  of  combustion,  J/kg. 

an  index  of  1  for  laminar  flow  and  of  0  for  turbulent  flow, 
ambient  gas  conductivity,  W/mK. 

distance  along  the  bed  from  the  origin  fixed  at  the  original 
ignition  front,  m. 
latent  heat  of  pyrolysis,  J/kg. 

pyrolyzate  transpiration  rate  per  unit  fuel  bed  width,  kg/ms 
pyrolyzate  mass  flux,  kg/m^s. 
an  index. 

upstream  heat  flow  rate  per  unit  fuel  bed  width,  W/m. 

combustive  heat  release  rate  per  unit  fuel  bed  width,  W/m. 

heat  flux,  W/m2, 

ambient  gas  temperature,  K. 

pyrolysis  temperature,  K. 

solid  surface  temperature,  K. 

flame  spread  velocity,  m/s. 

nondimensional  flame  spread  velocity. 

length  of  fuel  bed  between  pyrolysis  and  burnout  fronts,  m. 
location  of  pyrolysis  front,  m. 
ambient  gas  oxygen  mass  fraction, 
ambient  gas  kinematic  viscosity,  m?/s. 


<<>  a  nondimensional  ratio  of  heat  to  mass  fluxes  in  Sibuklin 

analysis. 

ps  solid  density,  kg/nv*. 

Subsection  8.4; 

A,B,C,D  Atallah's  forest  fire  spread  nondimensional  parameters. 
A,B,C,D  modified  At  a 11  a  h  parameters. 

Cp,  CpS  solid  fuel  specific  heat,  J/kg  K. 

d,  d0  fire  base  diameter,  m. 

E  plume  dynamics  nondimensional  parameter, 

f  Welker's  flame  shape  factor. 

F  fuel  bed  structure  parameter. 

Fr  Froude  number, 

g  gravitational  constant,  m/s2, 

g  Emmon's  firebreak  width,  m. 

Ti"  convective  heat  loss  coefficient,  W/m2K. 

hc  heat  of  combustion,  J/kg. 

H  Atallah's  flame  height,  m. 

H  Emmons/Shen  paper  strip  height,  m, 

k  entrai nment/mi xi ng  constant. 

K  a  constant  in  Welker's  work. 

Kg  gas  conductivity,  W/mK. 

K,  Ks  solid  conductivity,  W/mK. 

I  stick  height,  m. 

gas  phase  conduction  length,  m. 

Emmons/Shen  flame  height,  m. 


•j 

If 


fuel  loading  per  unit  area,  kg/m2. 
average  burning  rate,  kg/m^s. 
gaseous  fuel  feed  rate,  kg/s. 
an  index. 

fuel  moisture  content,  kg/kg. 

an  index  in  Vogel /Wi 1 1 iams  work,  1/4  for  laminar  flames  and 
0  for  turbulent  flames, 
an  index. 

Emmon's  burning  intensity  parameter. 

radiant  emission  intensity  per  unit  width,  W/m. 

ignition  energy,  J/kg. 

heat  output  rate  per  unit  length  of  a  line  fire,  W/m. 
heat  release  rate,  W. 

Vogel/Williams  downward  flame  propagation  speed  on  a  stick,  m/s. 
Steward's  stoichiometric  air/fuel  mass  ratio,  kg/kg. 

Emmon's  nondimensional  fire  spread  velocity. 

stick  spacing,  m. 

spacing  of  paper  strips,  m. 

time  to  Ignition,  s. 

ignition  temperature,  K. 

flame  temperature,  K. 

gas  temperature,  K. 

solid  initial  temperature,  K. 

pyrolysis  temperature,  K. 

ambient  gas  temperature,  K. 

Tf  -  Too,  K. 
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u  flame  spread  velocity,  m/s. 

Ug  fuel  feed  velocity,  m/s. 

utip  flame  gas  velocity  at  the  flame  tip,  m/s. 

V  nondimensional  flame  spread  velocity, 

v  wind  speed,  m/s. 

V  volumetric  fuel  feed  rate,  m3/s. 


w 

w 

wi 

x 

Xf 

y* 

z 

a 

as 

5 
3v 
e 

ef 

es 

0 

ua 

6 


p. 


Pb 

PO 


line  fire  width,  m. 
width  of  burning  zone,  m. 
ignition  strip  width,  m. 
distance  down  stream,  m. 
flame  tip  height,  m. 
flame  stand-off  distance,  m. 
height  from  floor,  m. 

Emmons'  fuel  bed  surface/volume  ratio,  1/m. 
solid  thermal  diffusivity,  ffl^/s. 

Rothermel's  porous  bed  packing  ratio, 
volumetric  expansion  coefficient,  1/K. 
an  effective  heating  number, 
flame  emissivity. 
solid  surface  emissivity. 
angle  of  tilt  of  a  wind-blown  flame, 
ambient  gas  dynamic  viscosity,  kg/ms. 
Rothermel's  propagation  flux  ratio. 
ps  solid  density,  kg/m3. 

Rothermel's  solid  bulk  density,  kg/m3, 
inlet  fuel  density,  kg/m3. 
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Pa*  p» 


ambient  air  density,  kg/m3, 
pf  flame  gas  density,  kg/m3. 

a  Rothermel's  fuel  particle  surface/volume  ratio,  1/m. 

a  Stefan-Bol tzmann  Constant,  W/m3k^. 

t  Emmons/Shen  paper  thickness,  m. 

t  forest  fuel  bed  height,  m. 

t  fuel  bed  thickness,  height  of  a  tree,  m. 

$  nondimensional  density  defect  parameter. 

4>s  slope  factor, 

ify,  wind  factor. 

Subsection  3.5: 

A  ceiling  +  Soffitt  skirt  area,  m?. 

Cpg  gas  specific  heat,  J/kg  K. 

Cpw  wall-lining  specific  heat,  J/kg  K. 

f  area  covered  by  a  flame  of  a  given  power,  m2/kW. 

h  convective  heat  transfer  coefficient  between  the  flame 

and  wal 1 ,  W/m2K. 

"h  convective  +  radiative  heat  transfer  coefficient  between 
the  hot  gases  and  wall-lining,  W/m3K. 

H0  window  height,  m. 

I  radiant  flux  imposed  in  a  heat  release  rate  calorimeter, 

k  a  window  flow  constant,  -  1.6  kg/m5/2s. 

Kw  wall-lining  conductivity,  W/mK. 

ma  air  flow  rate,  kg/s. 

mq  product  gas  outflow  rate,  kg/s. 


mv  fuel  feed  rate,  kg/s. 

q  measured  heat  release  rate,  W/m2. 

t  time,  s. 

Ta.T*  ambient  air  and  Initial  temperature,  K. 

Tg  gas  temperature,  K. 

Tw  wall  surface  temperature,  K. 

W0  window  width,  m. 

0  wall-lining  combustibility  parameter. 

y  ceiling  +  Soffitt  skirt  area  parameter. 

9g  nondimensional  gas  temperature. 

^  nondimensional  wall  surface  temperature, 

p*  wall-lining  density,  kg/m3. 

t  nondimensional  time. 

4>  nondimensional  fuel  feed  rate, 

t  nondimensional  air  flow  rate. 

4/  window  aspect  ratio, 

n  wall-lining  heating  character! sti c . 

Subsection  8,6: 

Ic  critical  irradiance  for  ignition,  W/m2. 

I0  emission  intensity  of  a  fire,  W/m2. 

Fn_f  radiation  view  factor  from  nth  window  flame  to  target  facade  f. 
n  window  number. 


N 


total  windows. 
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Section  9 
Subsection  9.1: 


A 

b 

c 

g 

Q 

r 

R 

Ro 

T 

To 

T„ 

T, 


source  energy  strength,  m^K/s. 
a  measure  of  plume  thickness,  m. 
plume  spreading  constant  =  dR/rlz. 
gravitational  acceleration,  m/s^. 
buoyancy  source,  m4/s3. 
radial  location,  m. 
radius  of  plume  at  any  height  z,  m. 
plume  source  radius,  m. 
temperature  at  any  (z,r),  K. 
source  temperature,  K. 
ambient  temperature,  K. 
maximum  temperature  at  any  z,  K. 
axial  velocity  at  any  (z,r),  m/s, 
source  axial  velocity,  m/s. 
maximum  axial  velocity  at  any  z,  m/s. 
nondimensi onal  axial  velocity, 
axial  location,  m. 

distance  of  virtual  point  source  below  the  finite  real 
source,  m. 

zceiling  ceiling  height  in  stable  atmosphere,  m. 
a  entrainment  constant. 

0  Rouse  constant  *  96/71. 

Y  stratification  number. 

9  nondimensi onal  temperature. 


'max 

w 

w0 

wmax 

W 

z 

zO 
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p  plume  gas  density,  kg/m3. 

pa,  ambient  gas  density,  kg/m3. 

Subsection  9.2: 

A  Morton's  nondimensional  radiation  free  path. 

C  Morton's  modified  Fronde  number. 

Cp  plume  fluid  specific  heat,  .3/kg  K. 

0  Morton's  ratio  of  radiative  loss  to  convective  transport. 

I  a  shape  factor. 

* 

nondimensional  ambient  pressure, 
q  power  of  the  buoyancy  source,  W. 

u  radial  velocity,  m/s. 

ue  entrainment  velocity,  m/s. 

c*,  weakly  buoyant  plume  entrainment  constant, 

a*  a  constant, 

Y  adiabatic  exponent. 

<  plume  gas  absorption  coefficient,  1/m. 

a  Stefan-Boltzmann  Constant,  W/m3K4. 

<t>  Murgai's  radiation  number. 

Subsection  9.3: 

L  fi re-base  scale,  mi 

r0  mean  radius  of  the  ring-fire-base,  m. 

U  nondimensional  radial  velocity. 

W  nondimensional  axial  velocity. 

<0  characteristic  velocity,  m/s. 
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Subsection  9.4; 


A  frontal  area  of  the  flying  brand,  m2. 

Cq  drag  coefficient, 

g  gravitational  constant,  m/s2, 

m  mass  of  the  firebrand,  kg. 
r  radius  of  the  firebrand,  m. 

r0  initial  radius  of  the  firebrand,  m. 

Re  Reynolds  number  of  flight, 

t  1 1  me ,  s . 

tb  burnout  time,  s. 

tg  time  of  fl ight,  s. 

Ui  a  nondimensi onal  parameter, 

wf  fall  velocity,  m/s. 

wf0  initial  fall  velocity,  m/s. 

8,  6  empirical  combustion  constants,  1/s2  and  m, 
C  density  ratio  p«/pwo* 

n  an  empirical  combustion  constant,  1/s2. 

XI, A2  length  ratios. 

ambient  air  density,  kg/m2, 
wood  density,  kg/m2, 
initial  wood  density,  kg/m2, 
nondimensional  time. 

Subsection  9.5: 

Cp  isobaric  heat  capacity,  J/kg  X. 

C v  isochoric  heat  capacity,  J/kg  K. 


Poo 

cw 

PWO 

X 


d  fire  diameter,  m. 

F  nondimensional  circulation  parameter, 

g  gravity,  m/s2, 

k  building  or  tree  height,  m. 

k*  surface  roughness  ratio  k/d. 

L  a  characteristic  length,  m. 

m  a  thermodynamic  parameter. 

P  pressure,  Pa. 

Fo  potential  pressure  of  the  atmosphere.  Pa. 

r  radius,  m. 

S  Ap/  Po 

t  t i me ,  s . 

v  velocity,  m/s. 

a  a  constant. 

e  nondimensional  gradient  of  potential  density  with  height 
p  density,  kg/m2. 

"p  potential  density,  kg/m2. 

"p0  potential  density  of  ambient  air,  kg/m2. 

Ap  po  “  p»  kg/m2, 

a  adiabatic  exponent,  Cp/Cv. 

$2  ambient  air  angular  velocity,  1/s. 

Subsection  9.6: 

C p co  ambient  air  specific  heat,  J/kg  K. 

Cpg  gas  specific  heat,  J/kg  K. 

d,dparticle,£  particle  diameter,  m. 


d0  particle  Initial  diameter,  m. 

dfoed*£  bed  diameter,  m. 

f  friction  factor, 

g  gravity,  m/s2. 

h  convective  heat  transfer  coefficient,  W/m2K. 

Hbed  bed  height,  m. 

j  an  index,  0,  1,  2  respectively  for  slab,  cylinder  and  sphere 

Kg  gas  conductivity,  W/mK. 

l  bed  dimension,  m;  also  particle  size,  m. 

i0  d0/2,  half-thickness  of  the  particle,  m. 

LSg  solid-to-gas  phase  change  latent  heat,  J/kg. 

m  mass  of  particle,  kg. 

m0  initial  mass,  kg. 

mp  burning  rate,  kg/s. 

rn^ax  maximum  burning  rate,  kg/s. 

m"  mass  flux,  kg/m2s. 

q“  heat  feed-back  flux,  W/m2. 

q"  radiant  heat  flux,  W/m2. 

rad 

Q  nondimensional  heat  flux, 

t  time,  s. 

t*  characteri sti c  pyrolysis  time,  s. 

Tb  burned  gas  temperature,  K. 

To,  hot  gas,  flame,  or  ambient  temperature,  K. 

Tp‘  pyrolysis  temperature,  K. 

T0  fuel  initial  temperature,  K. 


Tw  evaporating  fuel  surface  temperature,  K. 

^induced  Induced  wind  velocity,  m/s. 
oq  virgin  solid  fuel  thermal  diffusivlty,  m3/s 

6  fuel  bed  void  fraction. 

X  particle  shape  factor. 

p0  Initial  wood  density,  kg/m3. 

pc  final  char  density,  kg/m3. 

Pa,  ambient  air  density,  kg/m3. 
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